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A  HIGH  FREQUENCY  ANALYSIS  OF  ELECTROMAGNETIC  PLANE  WAVE 
SCATTERING  BY  A  FULLY  ILLUMINATED  PERFECTLY  CONDUCTING 

SEMI-INFINITE  CONE 


Keith  Dennis  Trott 
Captain,  U.S.A.F. 

Doctor  of  Philosophy  (Ph.D.) 

The  Ohio  State  University 
1986 

Professor  P.H.  Pathak,  Advisor 
127  pp. 

A  uniform  high  frequency  asymptotic  solution,  based  upon  the  physical  optics 
(PO)  approximation,  is  obtained  in  the  format  of  the  uniform  theory  of  diffraction 
(UTD)  to  describe  the  fields  diffracted  by  the  tip  of  a  semi-infinite,  perfectly  con¬ 
ducting  cone  when  it  is  fully  illuminated  by  an  electromagnetic  plane  wave.  The 
solution  is  expressed  in  terms  of  an  integral,  over  finite  limits,  which  can  be  nu¬ 
merically  integrated  without  difficulty.  The  solution  is  uniform  in  the  sense  that, 
unlike  the  previous  solution  based  upon  PO  which  is  valid  only  on  the  deep  shadow 
side  of  the  reflection  boundary,  and  becomes  singular  at  the  reflection  boundary, 
the  present  solution  remains  bounded  and  continuous  within  the  transition  region 
adjacent  to  the  reflection  boundary.  Outside  the  transition  region,  the  various  field 
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terms  in  the  integral  form  of  the  solution  can  be  interpreted  in  terms  of  the  geomet¬ 
rical  theory  of  diffraction  (GTD)  as  being  associated  with  the  total  field  consisting 
of  the  incident,  surface  reflected,  and  tip  diffracted  ray  fields.  It  is  noted  that  the 
PO  approximation  is  based  upon  the  use  of  the  geometrical  optics  (GO)  current  on 
the  cone.  A  correction  to  the  PO  tip  diffracted  field  is  investigated  based  on  the 
behavior  of  the  exact  eigenfunction  expansion  near  the  tip;  this  this  problem  needs 
further  study.  The  results  computed  from  just  the  uniform  asympototic  PO  solution 
compare  well  with  previously  published  results  for  narrow  angle  semi-infinite  cones.  <£~ — 

In  addition,  they  compare  well  with  measurement  as  well  as  with  an  independent 
moment  method  (MM)  solution,  for  the  scattering  by  a  finite  fiat-backed  cone.  In  , 

the  latter  case,  several  higher  order  wave  interactions  are  found  to  be  significant; 
one  such  is  an  interaction  between  the  tip  and  the  base  of  the  cone.  It  appears 
that  until  now,  there  were  no  expressions  in  the  open  literature  to  calculate  this 
tip-base  interaction  even  though  it  was  suspected  to  be  important.  The  present 
work  provides  the  expressions  necessary  to  calculate  this  interaction  and  also  to 
confirm  its  relative  importance.  _ _ 
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CHAPTER  I 
INTRODUCTION 

1.1  STATEMENT  OF  THE  PROBLEM 

A  uniform  high  frequency  asymptotic  solution,  based  upon  the  physical  optics 
(PO)  approximation  is  developed  to  describe  the  fields  diffracted  by  the  tip  of  a 
semi-infinite,  perfectly  conducting  cone  in  free  space  when  it  is  fully  illuminated  by 
an  electromagnetic  plane  wave.  This  solution  is  uniform  in  the  sense  that,  unlike 
the  exact  evaluation  of  the  far  zone  PO  integral  given  in  Chapter  II,  which  is  valid 
only  in  the  deep  shadow  side  of  the  reflection  boundary  and  becomes  singular  at 
the  reflection  boundary,  the  uniform  solution  presented  here,  remains  bounded  and 
continuous  within  the  transition  region  adjacent  to  the  reflection  boundary. 

This  solution  is  derived  using  asymptotic  methods  given  by  Felsen  and  Mar- 
cuvitz  [36].  The  form  of  the  solution  which  is  found  by  simply  applying  these 
techniques,  yields  a  uniform  asymptotic  theory  (UAT)  type  solution.  In  the  UAT 
form  [40],  the  total  field  is  written  as 

U  =  UG  +  U£  (1.1) 

in  which  the  Keller  diffracted  field,  u£,  is  based  on  Keller’s  geometrical  theory  of 
diffraction  (GTD)  where  one  writes  the  total  GTD  field  as 

U gtd  =  U?0  4-  Ujj,  (1.2) 

in  which  U9°  is  the  geometrical  optics  (GO)  field.  It  is  noted  that  becomes 
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singular  at  the  incident  shadow  boundary  (ISB)  and  the  reflection  shadow  bound¬ 
ary  (RSB).  The  U ®  in  (1.1)  is  a  modified  GO  field  that  eliminates  the  singularities 
of  (/£  at  the  ISB  and  the  RSB.  In  the  uniform  geometrical  theory  of  diffraction 
(UTD)  form  [40],  the  total  field  is  written  as 

Vutd  =  V9°  +  Vd  (1.3) 

where  XJd  is  a  uniform  diffracted  field.  It  is  noted  that  reduces  to  away  from 
the  shadow  boundary  transition  regions,  and  it  compensates  for  the  discontinuities 
of  U^0  at  the  shadow  boundaries.  Ur.Hke  UAT,  the  UTD  compensates  for  the 
singularities  of  t/jjf  at  the  ISB  and  RSB  in  a  multiplicative  manner  because  Ud  is 
expressed  as  a  product  of  and  a  transition  function  which  keeps  Ud  bounded  at 
the  ISB  and  RSB.  From  a  numerical  point  of  view,  the  latter  is  easier  to  implement 
in  computer  applications  than  the  additive  cancellation  of  the  singularities  in  the 
UAT  form.  Outside  the  transition  region,  both  the  UAT  and  UTD  forms  in  (1.1) 
and  (1.3)  reduce  to  Keller’s  GTD  form  in  (1.2). 

As  mentioned  previously,  the  uniform  solution  found  for  the  cone  starts  out 
in  a  UAT  form.  Steps  are  taken  to  rewrite  this  in  a  UTD  form.  As  a  first  step, 
a  quasi  UTD  form  is  found  for  the  scattered  field  that  is  expressed  in  terms  of  a 
relatively  simple  integral  over  finite  limits  which  can  be  integrated  numerically. 
The  integrand  is  written  as 

1  =  lref  +  hran  +  hiff  (14) 

where  the  integrands  denoted  by  I rej,  Itran,  and  ^diff  are  associated  with  the 
reflected  field,  the  transition  field,  and  the  Keller  type  tip  diffracted  field  respec¬ 
tively.  Far  from  the  transition  region,  Ifron  becomes  small,  and  the  remaining 
terms  in  (1.4)  can  be  interpreted  in  the  GTD  format  as  being  associated  with  the 
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surface  reflected,  and  tip  diffracted  ray  fields.  This  quasi  UTD  form  also  provides 
an  additive  correction  for  the  singularities  at  the  ISB  and  RSB  since  the  singularity 
in  I*//  *s  compensated  for  by  the  singularity  in  I*ran;  however,  it  is  desirable  to 
obtain  a  UTD  type  expression  for  the  tip  diffracted  field.  This  has  been  accom¬ 
plished.  In  the  development  of  the  latter  UTD  solution,  the  quasi  UTD  form  for 
the  scattered  field,  U,c,  is  rewritten  in  the  UTD  form  as 

Usc  =  Ur  +  Ud.  (1.5) 

where  U^0  is  the  reflected  component  of  the  GO  field,  and  is  the  uniform  tip 
diffracted  field.  This  form  gives  us  the  UTD  tip  diffracted  field  based  on  the  PO 
approximation.  This  UTD  tip  diffracted  field  is  valid  everywhere  away  from  the 
tip,  including  the  cone  surface.  It  should  be  noted  that  this  PO  based  tip  diffracted 
field  contains  a  factor  which  keeps  bounded  at  the  reflection  shadow  boundary, 
and  which  keeps  the  total  scattered  field  continuous  across  the  reflection  shadow 
boundary. 

The  present  uniform  solution  is  useful  for  analyzing  the  near  axial  radar  cross 
section  (RCS)  of  the  finite  flat-backed  cone.  In  the  RCS  analysis  of  the  finite  flat- 
backed  cone,  the  fields  diffracted  by  the  tip  travel  along  the  cone  surface  to  the 
base  and  diffract  from  the  edge  of  the  base  to  arrive  at  the  observation  point,  or 
travel  across  the  base  to  diffract  from  a  diametrically  opposite  edge  of  the  base  to 
arrive  at  the  observation  point,  and  vice  versa.  This  tip-base  interaction  requires 
one  to  calculate  the  tip  diffracted  field  on  the  illuminated  surface  of  the  cone. 
In  general,  the  illuminated  surface  lies  within  the  reflection  boundary  transition 
region;  therefore,  it  becomes  necessary  to  employ  a  uniform  solution  which,  as 
the  one  developed  here,  remains  valid  within  the  reflection  boundary  transition 
region.  It  appears  that  until  now,  there  were  no  expressions  in  the  open  literature 
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to  calculate  this  tip-base  interaction  even  though  it  was  suspected  by  [10]  to  be 


important.  The  present  work  provides  the  expressions  necessary  to  calculate  this 


interaction  and  also  to  confirm  its  relative  importance.  Indeed,  previous  work  on 


the  axial  RCS  of  a  flat-backed  cone  of  radius  a  for  the  range  of  8  <  ka  <  12,  has 


shown  an  almost  constant  3  dB  error  between  calculations  and  measurements  be¬ 


cause  the  tip-base  interaction  could  not  be  calculated.  That  error  is  now  removed 


by  the  present  work  thus  bringing  the  calculations  into  close  agreement  with  the 


measurements. 


Since  this  uniform  solution  allows  one  to  determine  the  value  of  the  field  on  the 


surface  of  the  cone  that  has  been  diffracted  by  the  tip;  therefore,  via  reciprocity, 


if  one  puts  an  infinitesimal  current  element  on  the  surface  of  the  cone,  this  same 


UTD  tip  diffracted  field  can  be  used  to  determine  the  contribution  to  the  radiated 


field  of  the  current  element  due  to  that  interaction.  An  example  of  how  reciprocity 


is  applied  is  given  in  Chapter  IV  in  terms  of  the  exact  eigenfunction  solution. 


It  is  noted  that  the  PO  approximation  is  based  upon  the  use  of  the  GO  current 


on  the  cone.  The  eigenfunction  expansion  is  investigated  and  the  dominant  term 


in  the  eigenfunction  expansion  is  isolated.  The  current  from  the  dominant  term 


exhibits  the  proper  singularity  at  the  tip.  It  is  compared  to  the  GO  current  for  the 


case  of  axial  incidence  and  it  is  also  shown  that  there  is  a  close  correlation  between 


these  two  expressions.  This  leads  one  to  believe  that  a  correction  to  the  PO  tip 


diffracted  field  could  be  constructed  using  an  expression  for  the  current  that  is 


the  matching  of  the  eigenfunction  current  near  the  tip  to  the  GO  current  far  from 


the  tip.  It  is  felt  that  this  could  provide  a  very  good  approximation  to  the  actual 


current  because  the  eigenfunction  expansion  yields  the  proper  tip  singularity,  and 


the  GO  current  dominates  far  from  the  tip  when  the  cone  is  fully  illuminated. 


1.2  MOTIVATION 


The  right  circular  cone  is  one  of  the  simplest  canonical  shapes  that  exhibits 
a  tip  discontinuity.  Little  work  has  been  done  regarding  tip  diffraction  in  the 
forward  scatter  region  toward  the  cone  surface.  This  information  is  needed  to 
accurately  study  the  scattering  from  the  cone  geometry,  and  thus  motivates  the 
present  research  effort. 

1.3  DISCUSSION  OF  PREVIOUS  RESEARCH 

There  has  been  a  wealth  of  work  conducted  on  the  cone  geometry;  both  finite 
and  semi-infinite.  The  following  summary  is  not  an  exhaustive  one.  Its  goal  is 
to  give  the  reader  an  indication  of  the  amount  of  interest  that  exists  in  the  cone 
geometry. 

Between  1952-1953,  Siegel,  et  al.  [1,2]  published  work  which  determined  scat¬ 
tering  from  a  cone  and  later  did  a  comparison  of  theory  to  experiment.  Their 
analyses  were  based  on  an  eigenfunction  expansion  that  is  slowly  convergent,  and 
special  summation  techniques  were  employed. 

In  1956,  Bailin  and  Silver  [3]  solved  the  boundary  value  problem  for  spheres 
and  cones.  They  also  used  the  radial  transmission  eigenfunction  expansion.  Pat¬ 
terns  were  plotted  for  circumferential  excitation  of  the  semi-infinite  cone.  A  few 
of  their  numerical  results  were  later  found  to  be  in  error  by  Pridmore-Brown  et 
al.  [35].  This  is  verified  by  the  present  work. 

In  1957,  Felsen  [4,5]  published  companion  articles  on  alternate  field  represen¬ 
tations  for  regions  bounded  by  spheres,  cones  and  planes.  In  these  papers,  he  uses 
the  angular  transmission  eigenfunction  expansion.  This  representation  converges 
faster  than  the  radial  transmission  version  used  by  Siegel,  et  al.  and  Bailin  and 
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Silver.  His  result  gave  the  tip  scattered  field  which  is  valid  only  on  the  shadow 
side  of  the  reflection  shadow  boundary  for  wide  and  narrow  angle  cones.  More¬ 
over,  his  eigenfunction  expansion  lends  itself  more  easily  to  an  asymptotic  high 
frequency  analysis.  In  1957,  Felsen  [6]  extracted  the  asymptotic  expansion  of  the 
tip  diffracted  wave  from  the  semi-infinite  cone  from  his  representation. This  result 
was  also  valid  only  on  the  shadow  side  of  the  reflection  shadow  boundary. 

In  1965,  Shultz,  et  al.  (7]  published  work  on  a  spherically  capped  finite  cone 
for  axial  incidence.  They  used  the  vector  wave  functions.  Their  results  are  valid 
up  to  ka=3.2,  where  a  is  the  radius  of  the  cone  base. 

Also  in  1965,  Bechtel  [8]  published  work  on  the  finite  cone.  Bechtel  used  the 
geometrical  theory  of  diffraction  (GTD).  The  upper  limit  in  his  study  was  a=2.87. 
He  neglected  tip  diffraction  and  did  not  account  for  the  tip-base  interaction. 

In  1969,  Bowman  et  al.  [33],  gave  a  complete  overview  of  most  of  the  work  done 
on  the  cone  geometry  up  to  that  point  in  time.  He  discussed  both  the  acoustic  and 
electromagnetic  cases.  This  treatise  was  used  to  compare  many  of  the  intermediate 
results  in  this  research. 

In  1971,  Senior  and  Uslenghi  [9]  published  work  on  the  backscatter  from  the 
finite  flat-backed  cone  using  GTD  and  keeping  the  first  two  terms  in  the  asymptotic 
expansion.  They  used  wedge  diffraction  coefficients  for  the  higher  order  diffractions 
across  the  base  of  the  flat-backed  cone.  There  is  an  error,  by  a  factor  of  2,  in  their 
treatment  of  grazing  incidence  which  was  later  corrected  in  1972  by  Burnside  and 
Peters  [10]  who  used  equivalent  currents  and  higher  order  diffractions  to  determine 
the  axial  radar  cross  section  (RCS)  of  the  finite  flat-backed  cone.  This  paper  also 
user  the  equivalent  wedge  diffraction  coefficients  for  higher  order  diffractions  across 
the  base  of  the  cone.  They  indicate  a  3  dB  discrepancy  between  measured  and 
calculated  data  for  a  cone  where  8  <  ka  <  12,  and  postulated  this  to  be  caused  by  a 


missing  diffraction  mechanism;  namely,  the  tip-base  interaction.  Their  conjecture 
is  verified  to  be  true  in  the  present  work. 

In  1977,  Kouyoumjian  [11]  published  work  on  the  high  frequency  backscatter 
from  fiat  backed  finite  cone.  This  paper  presents  both  the  acoustic  and  electromag¬ 
netic  solutions.  In  this  paper  the  first  attempt  is  made  to  account  for  the  tip-base 
interaction.  An  expression  is  used  for  the  tip  diffraction  that  is  not  strictly  valid 
beyond  the  reflection  shadow  boundary;  however,  the  results  are  reasonably  ac¬ 
curate  for  large  ka  but  there  is  still  more  than  a  one  dB  difference  for  smaller  ka 
values. 

From  1977-1979,  Felsen  and  Chan  published  a  series  of  papers  [12,13,14]  in 
which  the  integral  form  of  the  Green’s  function  for  the  cone  was  determined.  The 
first  of  these  extracted  the  creeping  wave  from  two  points  located  on  the  surface 
far  from  the  tip.  In  the  latter  two  papers,  their  Green’s  function  was  transformed 
to  the  time  domain  to  investigate  transient  behavior. 

In  1984,  Wang  and  Mitschang  [15]  published  a  hybrid  moment  method  (MM) 
and  physical  theory  of  diffraction  (PTD)  solution  for  the  finite  circular  and  el¬ 
liptical  cones.  They  commented  in  their  study  that  the  tip-base  interaction  was 
negligible;  however,  that  comment  is  shown  to  be  untrue  in  the  present  work. 

It  has  come  to  our  attention  that  Molinet  [39]  has  also  independently  worked 
on  the  problem  of  scattering  from  the  cone  tip  using  PO  and  high  frequency  tech¬ 
niques  from  Felsen  and  Marcuvitz  [36].  He  has  treated  the  special  case  of  axial 
incidence  [39].  His  results,  which  were  determined  independently,  were  used  to 
verify  the  general  solution  given  here;  however,  unlike  the  present  UTD  form  of 
the  solution,  he  keeps  his  high  frequency  solution,  also  based  on  PO,  in  the  form 
of  uniform  asymptotic  theory  (UAT). 

This  work  is  devoted  to  determining  the  value  of  the  wave  that  is  excited  at 


the  tip  via  diffraction  and  travels  along  the  cone  surface  until  diffracted  by  the 
edge  at  the  base  of  a  finite  cone.  A  uniform  high  frequency  asymptotic  solution, 
based  upon  the  physical  optics  (PO)  approximation,  is  obtained  to  describe  the 
fields  diffracted  by  the  tip  of  a  semi-infinite,  perfectly  conducting  cone  when  it 
is  fully  illuminated  by  an  electromagnetic  plane  wave.  The  solution  is  uniform 
in  the  sense  that,  unlike  the  far  field  solution  based  upon  PO  which  is  valid  only 
in  the  deep  shadow  side  of  the  reflection  boundary,  the  present  solution  can  be 
employed  even  within  the  transition  region  adjacent  to  the  reflection  boundary  up 
to  the  cone  surface.  The  results  computed  from  the  present  uniform  asymptotic  PO 
solution  are  shown  to  compare  well  with  previously  published  results  for  narrow 
angle  semi-infinite  cones  [5].  In  addition,  the  results  obtained  for  the  RCS  of  the 
finite  flat-backed  cone  are  found  to  compare  closely  with  measurements  and  with 
an  independent  moment  method  (MM)  solution.  It  should  be  noted  that  the  UTD 
tip  diffracted  field  on  the  cone  surface,  found  in  the  present  work,  can  be  used  in 
conjunction  with  reciprocity  to  determine  the  radiation  from  a  source  on  the  cone 
surface. 

1.4  FORMAT  OF  THE  DISSERTATION 

The  format  for  this  dissertation  is  as  follows.  In  Chapter  II,  an  exact  evalua¬ 
tion  of  the  far  zone  field  scattered  from  the  cone  is  developed  using  Physical  Optics 
(PO).  If  this  result  were  generalized  to  the  near  zone,  it  would  be  a  non-uniform 
expression  since  the  expressions  for  fields  are  singular  at  the  reflection  shadow 
boundary.  Furthermore,  it  is  not  valid  on  the  lit  side  of  the  shadow  boundary. 
The  reasons  for  this  restriction  will  be  discussed.  This  is  done  for  plane  wave  inci¬ 
dence  in  the  special  case  when  the  cone  is  fully  illuminated.  This  field  expression 
is  then  used  to  determine  the  scattering  cross-section,  a.  It  is  shown  that  this  field 


expression  yields  the  familiar  radar  cross  section  from  PO,  opo,  as  shown  in  [33]. 
Moreover,  the  field  expression  reduces  to  that  given  by  Goryainov  [37]  for  axial 
incidence.  It  appears  that  an  explicit  field  expression  for  the  case  of  non-axial 
incidence,  but  with  full  illumination,  has  not.  been  published  previously  although 
the  Op0  has  been  obtained  previously. 

In  Chapter  III,  a  uniform  solution  based  upon  PO  is  developed  for  scattering 
from  a  fully  illuminated  semi-infinite  cone.  In  the  PO  radiation  integral  over 
the  cone  surface,  the  free-space  Green’s  function  with  an  exact  kernel  is  used.  The 
integration  is  performed  asymptotically  in  the  radial  coordinate  (r')  along  the  cone 
generator  using  a  uniform  saddle  point  method  which  takes  into  consideration  the 
effect  of  the  saddle  point  near  the  end  point  (tip).  The  remaining  azimuthal  (^') 
integration  is  then  done  numerically.  It  is  found  that  the  saddle  point  migrates  to 
the  end  point  (tip)  as  a  function  of  <J>'. 

This  answer  is  first  written  in  a  quasi  uniform  theory  of  diffraction  (UTD) 
form.  It  includes  a  surface  reflected  field  term,  a  tip  diffracted  field  term,  and  a 
transition  term  which  vanishes  in  the  deep  shadow  region.  The  transition  term 
keeps  the  total  solution  bounded  and  continous  across  the  reflection  boundary. 
In  general,  the  transition  region  extends  throughout  the  lit  region  for  the  fully 
illuminated  cone.  This  quasi  UTD  solution  was  subsequently  rewritten  in  a  UTD 
form  such  that  a  UTD  type  tip  diffracted  field  has  been  found  using  the  PO 
approximation.  One  of  the  goals  of  this  research  is  to  develop  an  expression  that 
can  be  used  to  calculate  the  contribution  by  the  wave  diffracted  from  the  tip.  In 
the  case  of  finite  flat-backed  cone,  the  tip  diffracted  wave  travels  along  the  cone 
and  can  diffract  from  the  edge  at  the  base  of  the  cone.  The  UTD  tip  diffracted  field 
has  been  determined  for  the  semi-infinite  cone,  but  since  high  frequency  diffraction 
is  a  local  phenomena,  this  tip  diffraction  is  the  same  for  either  the  finite  or  the 


semi-infinite  cone. 

To  date  there  has  been  an  almost  constant  3  dB  difference  between  analysis 
and  experiment  for  axial  backscatter  from  a  moderately  large  (8  <  ka  <  12) 
flat-backed  finite  right  circular  cone  in  which  ”a”  is  the  radius  of  the  base  and 
k  is  the  free  space  wave  number.  It  is  shown  that  including  the  aforementioned 
tip-base  interaction  improves  the  results  and  brings  the  analysis  in  line  with  the 
experiments. 

In  Chapter  IV,  an  exact  eigenfunction  expansion  is  derived  for  the  fields  using 
the  vector  wave  functions  M  and  N  described  in  [21,32].  From  this,  expressions  are 
determined  for  the  surface  current  induced  on  the  cone.  When  the  small  argument 
form  of  the  Bessel  function  is  applied  to  the  series,  it  can  be  shown  that  there 
is  a  dominant  term  in  the  tip  region.  This  term  represents  the  tip  discontinuity 
and  models  the  true  behavior  of  the  current  near  the  tip.  For  the  case  of  axial 
incidence,  the  actual  tip  current  and  the  GO  based  current  are  compared  and  it  is 
shown  that  the  two  currents  are  very  similar.  This  leads  one  to  believe  that  these 
currents  could  be  combined  to  yield  a  more  accurate  approximation  to  the  current 
which  could  then  be  used  to  determine  the  tip  diffracted  field  using  the  methods 
in  Chapter  III. 

In  addition,  the  vector  wave  functions  are  employed  to  compute  the  radiated 
fields  from  a  circumferential  slot  on  the  cone.  Baiiin  and  Silver  [3]  have  published 
patterns  for  several  configurations.  Pridmore-Brown  et  al.  [35]  have  also  published 
patterns  for  circumferential  slots.  The  Pridmore-Brown  paper  claims  that  there 
are  several  discrepancies  in  the  results  of  Baiiin  and  Silver.  To  resolve  this  issue, 
patterns  are  plotted  based  on  the  present  work  using  the  same  parameters.  It 
should  be  noted  that  the  present  expressions  using  the  vector  wave  functions  are 
identical  to  those  given  in  [3],  although  derived  using  a  totally  different  approach. 


Moreover,  the  present  expressions  appear  vastly  different  from  those  in  [35].  The 
case  in  point  is  for  ka=10,  and  9q  =  170°.  The  present  plots  agree  with  those  in 
[35]  and  are  very  different  from  those  in  [3].  It  should  be  noted  that  the  present 
plots  do  agree  with  that  in  [3]  for  ka=7r,  and  Oq  =  165°.  One  could  attribute  the 
discrepancy  to  the  level  of  sophistication  of  computers  then  and  now. 

Appendix  A  gives  the  expressions  used  in  the  uniform  saddle  point  approx¬ 
imation  derived  in  Chapter  III.  The  derivation  of  the  vector  wave  functions  and 
the  dyadic  Green’s  function  used  in  the  present  work  are  detailed  in  Appendices 
B  and  C.  An  important  aspect  of  this  problem  is  the  heavy  dependence  on  special 
functions.  There  is  also  the  requirement  to  compute  the  incomplete  normalization 
integrals.  Even  more  critical  is  the  determination  of  the  eigenvalues  associated 
with  the  cone  problem.  This  requires  the  eigenvalues  that  are  solutions  to 

P™(cos0o)  =  0  and  (1.6) 

^P„m(cosl?o)  =  0.  (1.7) 

For  the  cone  problem,  the  eigenvalues  u  and  7  are  non-integer;  therefore,  a  special 
effort  was  required  to  derive  expressions  that  could  be  programmed  to  yield  results 
for  any  angle,  order  and  degree.  Programs  were  written  to  calculate  these  and  were 
compared  to  published  results  in  [16].  Integral  forms  for  the  associated  Legendre 
polynomials  found  in  [17]  were  used  for  these  functions.  These  integral  forms  were 
also  used  to  determine  the  eigenvalues.  These  eigenvalues  compare  very  well  with 
those  published  in  [2,31].  Appendix  D  discusses  the  derivations  of  the  associated 
Legendre  polynomials  used  in  this  dissertation. 

The  e]wt  time  convention  is  assumed  and  suppressed  throughout  the  following 
development. 
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CHAPTER  II 


EXACT  EVALUATION  OF  FAR  ZONE  FIELDS  USING  PHYSICAL 

OPTICS 

2.1  INTRODUCTION 

In  this  chapter,  the  field  scattered  by  the  cone  using  the  Physical  Optics  (PO) 
approximation  is  determined.  The  PO  radiation  integral  which  is  valid  in  the  far 
zone  of  the  tip  is  used  and  the  discussion  is  specialized  to  the  case  when  the  cone 
is  fully  illuminated.  This  is  true  when  0  <  it  —  Oq  =  a,  where  the  angles  are  shown 
in  Figure  1. 

It  will  be  shown  that  the  expressions  for  the  scattered  field  found  in  this  chap¬ 
ter  are  not  valid  at  or  beyond  the  reflection  shadow  boundary  due  to  restrictions 
placed  upon  the  integration  techniques.  The  integral  is  over  the  cone  surface  which 
is  described  in  terms  of  the  radial  (r1)  and  azimuthal  [<f>')  coordinates.  The  inte¬ 
gration  in  (r')  is  over  infinite  limits;  therefore,  these  restrictions  are  placed  upon 
the  integral  to  insure  convergence.  The  integration  has  been  performed  two  ways. 
The  first  approach  is  shown  in  detail  within  this  chapter.  It  was  also  done  in  the 
reverse  order  holding  ( <f> ')  fixed,  integrating  over  (r'),  then  performing  the  integra¬ 
tion  in  (<f> )  using  contour  integration.  Bo'.h  techniques  yield  the  same  answers  and 
the  same  restrictions. 

Since  these  expressions  for  the  scattered  field  are  valid  only  on  the  shadow 
side  of  the  reflection  boundary,  they  represent  the  fields  diffracted  by  the  tip.  If 


Figure  1:  Geometry  for  Plane  Wave  Incidence 

this  were  generalized  to  the  near  zone  case,  then  this  solution  would  be  considered 
non-uniform  since  it  is  singular  at  the  reflection  shadow  boundary. 

2.2  CALCULATION  OF  THE  SCATTERED  FIELD 


Figure  1  shows  the  geometry  for  a  plane  wave  incident  on  the  cone.  The 
electric  field  associated  with  the  incident  plane  wave  at  the  conical  surface  is  given 


EJ(r')  =  j  j  e>**-'(sin0O8i 


sin  3  cos  i^+cos  6q  cos  3) 


where  the  bracketed  term  represents  the  two  possible  polarizations.  The  variables 
(r')  and  (<£')  represent  the  radial  and  azimuthal  variation  of  the  source  coordinates 
as  shown  in  Figure  1.  It  is  assumed  that  the  plane  of  incidence  is  the  x-z  plane. 
For  this  configuration,  the  incident  magnetic  field  H5  field  is  given  by 


H*(r')  =  |  |  ybejfcr'<8in0o 


sin  3  cos  4>'+ cos  8q  cos  3) 


In  the  PO  approximation,  the  radiation  integral  for  the  scattered  field  is  eval¬ 
uated  after  employing  the  geometrical  optics  (GO)  approximation  for  the  surface 


current.  This  GO  current  is  given  by 

j  (r')  =  /  2 n'  x  H‘(r')  in  the  lit  region,  and  (2.3) 

8  v  '  (  0  in  the  shadow  region.  ' 

However,  the  entire  cone  is  illuminated;  therefore,  there  is  no  shadow  region  on 
the  surface,  and  the  integration  is  performed  over  the  entire  cone  surface. 

The  outward  normal  to  the  cone  is  given  by 

n'  =  —O'  =  -icos0ocos<y  —  y  cos  sin  +  2sin0o-  (2-4) 

■  9'=(>o 

Therefore,  the  expression  for  the  scattered  electric  field  in  the  far  zone  of  the  tip 


is  given  by 


e-)kr- 

E8(r)  ~  jk - f  x  f  x 


fU  f^-ejkr'qtt'K,s\n  90dr'd^  (2.5) 

Jo  Jo  2tt 


where 


isin^o  +  2  cos  9q  cos  $  for  9'  incidence,  and 

f(^>  )  —  '  x  cos  0q  sin /?  sin  +  2  cos  0q  cos /?  sin  <£'  for  <j>'  incidence 
+y  (sin  0O  cos  0  -  cos  6 q  sin  0  cos  <f>') 


q{4> )  =  sin  0q  [(sin  0cos  </>  +  sin /?)  cos(// +  sin  0sin  ^sini//] 


+  cos  0O  [cos  9  4-  cos  0] . 


where 


One  can  rewrite  this  as  follows: 


q[4>)  =  acos(<f>'  —  6)  +  7 


-7  =  cos  0o  [cos 0  +  cos/?) 
a  cos  6  =  sin  9q  (sin  9  cos  </>  +  sin  0)  and 


(2.10) 
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a  sin  6  =  sin  s^n  ^  s*n 


(2.11) 


From  Equations  (2.10)  and  (2.11),  it  is  found  that 

a  =  sin  0O  \J sin2  0  +  sin2  0  +  2  sin  6  sin  0  cos  <f>. 
The  expression  for  the  scattered  field  becomes 

e-jkr 

EB(r)  =£  jk - sin^o^  *  r  * 

*  2tt  t(A> 


(2.12) 


r'e]hr'  f  *  fi^le]kar'coB^'-S)d(i)ldrl  (2 .13) 

Jo  Jo  27r 


Concentrating  on  the  <f>*  integration,  the  following  results  [21]  are  used: 

— ~ —  f  cos  n4>e]XCOi^~^  )d<f>'  =  cos  n<f>Jn(x) 

2 *]n  Jo  V  ' 


(2.14) 


and 


f  sin  n<f>  eixco*^  ^  ^d(f>  =  sin  n<})Jn{x).  (2-15) 

1*3  Jo 


For  either  polarization,  there  will  be  linear  combinations  of  the  same  three  basic 
integrals.  Using  the  above  results,  they  are  given  by 

*2n 
1 0 


—  f  *  Jkar'cos {^-6)^1  =  JQ(kar>)  (2.16) 

2?r  J o 

—  f  cos  <f>' ejkar  =  j  cos  6  Ji{kar')  (2.17) 

2? r  J o 


and 


—  [  sin  <t>' ejkaT  C08^  ^ =  j  sin  6  Ji(kar').  (2.18) 

2n  Jq 

Once  the  <f>'  integration  has  been  performed,  what  remains  are  integrals  of  the 


form  : 


and 


ir' = 


(2.19) 


(2.20) 
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The  answers  to  these  were  found  by  applying  the  following  results  [38]: 


r  -°*r  la,  + 1/2) 

Jo  '  J*'(/,x)x  -  vSF(ffl2  +  ^,»+./J 


(2.21) 


«  axJv{0x)x 


„+l  _  go(2gn>  j-  3/2) 


(2.22) 


/o  '  ^F(tt2  +  ^S/2 

where  i?e(a)  >  |/m(/?)|.  It  has  been  assumed  that  k  is  slightly  lossy  so  this 


condition  yields  —7  >  a.  Substituting  =  ir  —  a  into  the  expressions  for  7  and  a 


yields  the  following  restriction  on  the  integral: 


cos  a  [cos  9  +  cos  0}  >  sin  ay  sin2  9  ■+■  sin2/?  +  2  sin  9  sin/?  cos  <f>.  (2.23) 


This  expression  will  be  discussed  later  after  the  shadow  boundaries  have  been 


discussed. 


Defining  Eg  as  the  9  component  of  the  field  generated  by  a  O’  source  and  so 
forth,  then  using  the  following  identities 


f  x  f  x  x  ~  —9  cos  9  cos  (f>  +  ^  sin  <f> 


(2-2-1) 


f  x  f  x  y  ~  —  9cos9s\n<f>  —  <t>cos(j> 


(2.25) 


f  x  r  x  z  =  9s\n9, 


(2.26) 


then  simplifying  the  expressions  yields: 


—jlc  r  2 

Eg  =  J~j~ — [[l  +  cos0cos/?]cos<£  +  sin0sin/?] 


(2.27) 


ai  ,e  •;fcrtan2a  .  .  , 

K  =  -J-r~  — J7^-[cos0  +  cos/?]sin^ 


^  J  kr  n3/2  1  ^  ^  j  —  ^ 

2?^  =  [[l  "1"  cos  9  cos  /?]  cos  <f>  +  sin  9  sin /?] 


(2.28) 


(2.29) 


a'  ,e  jfcrtan2ar  „  .  . 

E%  =  [cos  g  +  cos  0}  sin  <t> 


(2.30) 


where  II  is  given  by 


II  =  [cos0  +  cos/?}2  —  tan2 a  ^sin20  -f  sin2/?  +  2  sin  0  sin/?  cos  .  (2.31) 

It  should  be  noted  that  tan2  0q  has  been  replaced  with  tan2  a  in  all  the  expressions 
above  since  a  =  n  —  0q. 

For  the  sake  of  completeness,  it  is  also  noted  that  the  scattered  magnetic  field 
in  the  far  zone  of  the  tip  is  given  by 


HB  ~  r  x  K0E8. 


(2.32) 


The  field  expressions  will  have  a  singularity  when  the  denominator  goes  to 
zero.  These  singularities  occur  at  the  shadow  boundaries.  The  reflection  shadow 
boundary  that  corresponds  to  the  specular  line  of  scattering  from  the  entire  cone 
generator  occurs  in  the  incident  plane,  i.e.,  <f>  =  0.  Setting  <t>  =  0,  II  can  be  written 
as 

4  cos2(^  -— )  cos2(— — — — )  —  tan2  ccsin2(— ^— )  .  (2.33) 

Setting  this  equal  to  zero  yields  a  condition  on  0  +  0  given  by 


+  0  —  ■k  -  2a  —  20q  ~  7r- 


(2.34) 


This  is  the  reflection  shadow  boundary  (RSB).  The  RSB  also  occurs  elsewhere, 
and  its  position  varies  as  a  function  of  <f>. 

Setting  <f>  =  7r  yields  the  incident  shadow  boundary  (ISB)  given  by 


9  +  0  =  7T. 


(2.35) 


Going  back  to  equation  (2.23),  the  restriction  on  the  integral,  and  setting 


<t>  =  0,  yields 


/  9  +  /?. 

cos(a  H - )  >  0 

v  2  ’ 


(2.36) 


from  which  it  is  noted 


6  +  0  <  n  -  2a.  (2.37) 

Setting  <f>  =  7T,  yields 

0  +  0  <  7T.  (2.38) 

One  should  recognize  these  values  as  the  RSB  and  ISB,  respectively.  The  incident 
shadow  boundary  poses  no  problem  in  this  case.  However,  there  is  a  phase  shift 
passing  through  the  RSB  and  the  expressions  are  not  valid  beyond  this  boundary. 
The  expressions  will  be  plotted  beyond  the  RSB  by  introducing  the  phase  shift  that 
is  caused  by  passing  through  this  boundary.  The  plots  indicate  that  the  behavior 
is  reasonable  for  axial  incidence;  however,  it  shows  rapid  flueuations  on  the  lit  side 
of  the  reflection  shadow  boundary  for  non-axial  incidence. 

As  mentioned  previously,  these  fields  are  only  valid  in  the  shadow  region. 
They  do  not  account  for  the  surface  reflected  field  or  the  transition  field;  therefore, 
the  scattered  field  given  is  the  tip  diffracted  field. 

It  is  of  interest  to  compute  the  scattering  cross-section  from  these  expressions. 
The  scattering  cross-section  is  defined  by 

|E8|2 


a  —  lim  4n r2- — t-t  =  lim  4n r2  [|Ei|2  +  \Ei\2 
r— oo  EM2  r— oo  l1  e'  1  <t>' 


(2.39) 


|ET 

In  either  case,  after  simplification 

Opo  = - “Yjj —  [ [  1  •+-  cos  8  cos  /?]  -f-  cos  (f>  sin  0  sin  /?]2  (2.40) 

where  II  is  given  by  Equation  (2.31).  This  expression  for  apo  is  identical  to  that 
given  previously  by  [33].  Moreover,  Goryainov  has  done  this  for  axial  incidence. 
His  expressions  are  given  in  [37]  by: 


m  ,e~]kr 

Eg  -  j— — cos  <f>L(9,  a)  and 
E$  =  sintf>L(0,a) 


(2.41) 

(2.42) 
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where 


,  sin  acosft 

=  ,  3/2 ' 

4  cos  j  (cos(a  +  |)  cos  (a  -  j)  j 

(2.43) 

The  exact  PO  expressions  starting  with  Equation  (2.27),  for  0'  incidence  could  be 

written  as: 

/  jfcr 

Eg  =  j  ^  cos<t>Tl(0,a,0)  +  T2(0,  a,0)  and 

(2.44) 

El  =  J  kr  sin<j>T3(0,a,0) 

(2.45) 

where 

♦  /y 

Tl(0,a,0)  =  ^3/2  (jl +cos0cos/?]j 

(2.46) 

7. 

♦  on“  yy 

T2(0,a,0)  =  ^  [sinflsin/?]  and 

(2.47) 

T3(0,a,0)  =  [cos 0  + cos/?] 

(2.48) 

Evaluating  the  expressions  for  0  —  0,  yields 

T2(9,  a,0)  =  0  (2.49) 

and 

Tl(0,a,O)  =  T3(0,  a,0)  =  L(0,a).  (2.50) 

That  is,  the  expressions  reduce  to  that  given  by  Goryainov  for  axial  incidence. 

In  some  applications,  the  individual  scattered  field  expressions  are  more  useful 


than  the  composite  apo  expression.  It  is  important  to  note  that  the  scattered  field 
expressions  and  hence  op0  satisfy  reciprocity. 


CHAPTER  III 


A  UNIFORM  ASYMPTOTIC  ANALYSIS  OF  CONE  TIP 
DIFFRACTION  BASED  ON  PHYSICAL  OPTICS 

3.1  INTRODUCTION 

The  exact  far  zone  PO  integral  calulated  in  the  preceeding  chapter  was  sin¬ 
gular  at  reflection  shadow  boundary  (RSB)  and  not  valid  on  the  lit  side  of  the 
reflection  shadow  boundary.  In  this  chapter,  the  solution  for  the  scattering  by 
a  fully  illuminated  cone  is  improved  by  making  it  uniform  across  the  reflection 
shadow  boundary.  To  this  end,  an  exact  kernel  for  the  Green’s  function  is  used  in 
the  PO  radiation  integral  for  the  scattered  field  in  which  the  integration  is  taken 
over  the  GO  currents  on  the  cone  surface.  This  integral,  which  is  valid  in  the  near 
zone  of  the  tip  but  not  too  close  to  the  tip,  will  be  approximated  asymptotically 
in  the  radial  coordinate  (rf)  along  the  generator  of  the  cone  using  the  saddle  point 
method.  The  saddle  point  can  migrate  to  the  end  point  of  integration  at  the  tip  of 
the  cone;  therefore,  special  techniques  are  used  to  account  for  this.  The  integration 
in  the  remaining  azimuthal  coordinate  is  done  numerically.  This  is  done  only 
for  O'  incidence.  The  4>'  case  would  follow  in  a  similar  manner. 

To  determine  the  field  on  the  cone  surface,  a  uniform  solution  is  required. 
There  are  two  common  types  of  uniform  solutions.  These  will  be  discussed  briefly. 
As  discussed  in  the  Introduction,  the  total  field  in  the  UAT  form  [40]  is  written  as 

U  =  UG  +  U]J  (3.1) 
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in  which  the  Keller  diffracted  field,  u£,  is  based  on  Keller’s  geometrical  theory  of 
diffraction  (GTD)  where  one  writes  the  total  GTD  field  as 


vgtd  =  uSo  +  u  ^  (3.2) 

in  which  U^0  is  the  geometrical  optics  (GO)  field.  It  is  noted  that  Uj£  becomes 
singular  at  the  incident  shadow  boundary  (ISB)  and  the  reflection  shadow  bound¬ 
ary  (RSB).  The  UG  in  (1.1)  is  a  modified  GO  field  that  eliminates  the  singularities 
of  u£  at  the  ISB  and  the  RSB.  In  the  uniform  geometrical  theory  of  diffraction 
(UTD)  form  [40],  the  total  field  is  written  as 

vutd  =  +  u*  (3>3) 

where  is  a  uniform  diffracted  field.  It  is  noted  that  Ud  reduces  to  away  from 
the  shadow  boundary  transition  regions,  and  it  compensates  for  the  discontinuities 
of  JJgo  at  the  shadow  boundaries.  Unlike  UAT,  the  UTD  compensates  for  the 
singularities  of  U£ f  at  the  ISB  and  RSB  in  a  multiplicative  manner  because  is 
expressed  as  a  product  of  Uf  and  a  transition  function  which  keeps  bounded  at 
the  ISB  and  RSB.  From  a  numerical  point  of  view,  the  latter  is  easier  to  implement 
in  computer  applications  than  the  additive  cancellation  of  the  singularities  in  the 
UAT  form.  Outside  the  transition  region,  both  the  UAT  and  UTD  forms  in  (3.1) 
and  (3.3)  reduce  to  Keller’s  GTD  form  in  (3.2). 

The  solution  developed  here  is  in  terms  of  a  finite  integral.  The  integrands  for 
this  integral  are  first  developed  into  a  quasi  UTD  form;  however,  a  UTD  expression 
for  the  tip  diffracted  field  is  desired.  To  this  end,  knowledge  of  the  behavior  in  the 
vicinity  of  the  reflection  shadow  boundary  is  utilized  and  the  quasi  UTD  form  is 
rewritten  in  a  UTD  form  that  isolates  the  UTD  tip  diffracted  field. 


(r  x 


Figure  2:  Geometry  for  Plane  Wave  Incidence 

3.2  FORMULATION  OF  THE  NEAR  FIELD  PO  RADIATION  IN¬ 
TEGRAL 

The  expression  for  the  GO  current  for  the  case  being  considered  is  given  by 
Jgo^,^)  =  2Yoe]kr'cosl(xsin  a.  —  zcosacos<j>')  (3.4) 


where 


cos  7  =  sin  a  sin  0  cos  <t>  —  cos  a  cos  0. 


The  magnetic  vector  potential  A  defined  by  HB  =  V  x  A  associated  with  this  Jgo 


is  given  by 


r  2n  r  oo  t~JKn 

A(r )  =  JQ  d<t>  Jq  Jgo(»■^  <f>')  sin  adr' 


e-jkR 

i  j.i\e  J. 


R  =  |r  -  r'|  =  \Jr2  +  r2  -  2 rr'  cos  <; 


cos  f  =  sin  a  sin  0  cos (<j>  -  4?)  ~  cos  Q  cos  ® 


Using  the  fact  that  H®  =  VxA  the  scattered  magnetic  field  is  given  by 

r  2n  roo  e~jkR 

H"(r)  =  V  x  /  d<j>'  Jgo(r', <*')^rr' 8in <*dr'  (3.9) 

JO  Jo 

Since  the  curl  is  operating  on  the  unprimed  or  observation  coordinates,  one  can 
move  this  operation  under  the  integral.  Only  the  Green’s  function  is  a  function  of 
r;  therefore,  after  applying  a  vector  identity, 

V  x  (jgo{r',d>')G0(R))  =VG0(R)  X  Jgo  ('V)  (3.10) 

where 

e—jkR 

C°'R>  “  1ST 

and 

d  1  e~]kR  - 

VGoW  =  =  -(/*  +  (3.12) 

since 

VR  =  R  (3.13) 


(3.11) 
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(3.14) 


Substituting  into  (3.9),  one  obtains 

HB(r)  =  [2*  d<j>'  /"0Of(r^^,)e^r'-^dr, 

Jo  Jo 

where 

f(r',<t>')  =  -r'  (jk  -f  i).R  x  (xsina  -  z  cos  a  cos  <i>)  (3.15) 

2  7r  R  R 

and 

q(r' ,  <t>)  =  r*  cos-y  -  \J r2  +  r*2  -  2rr'  cos  f .  (3.16) 

Concentrating  on  the  inner  integral  on  r*  in  Equation  (3.9),  one  finds  that 

1(1^,^)=  /°0f(r',^)^r''^dr'.  (3.17) 

Jo 

The  function  f  is  slowly  varying  in  the  interval  and  the  exponential  function  is 
oscillatory  with  a  stationary  phase  point.  It  is  therefore  clear  that  the  dominant 
contribution  to  the  integral  in  (3.17)  comes  from  the  stationary  point  at  r'  =  ra 
and  from  the  end  point  or  cone  tip  at  r'  —  0.  Rather  than  use  the  stationary  phase 
method,  it  is  more  convenient  to  employ  the  saddle  point  technique  to  evaluate 
not  only  the  leading  terms  but  also  the  higher  order  terms  corresponding  to  the 
end  point  and  saddle  (stationary)  contributions  which  are  important  in  this  case. 
It  is  noted  that  it  is  usually  more  difficult  to  obtain  higher  order  terms  from 
the  stationary  phase  approximation.  The  saddle  point  method  requires  one  to 
integrate  along  the  steepest  descent  path  through  the  saddle  point.  Transforming 
to  the  steepest  descent  path  as  in  [36]  requires  the  following  definitions: 

<l(r',<t>')=q(rs,4>')-s2  (3.18) 

and 

G(s)=f(r',*')^  (3.19) 


where,  in  addition,  one  defines 


sa  =  ±\y/q{ra,4t)  -  q{0,<t>')\  rs  >  0. 


(3.20) 


G  has  not  been  written  as  a  vector  for  convenience  even  though  it  is  defined  in 
terms  of  a  vector  function  in  that  it  is  not  relevant  to  the  derivation. 


Next,  G(s)  is  written  as  in  [36]  as 


G(s)  =  G(0)  +  (G(s)  -  G(0)). 


(3.21) 


So  the  integral  I  becomes 


_  Jkq(rt,4>‘ 


I  =  e 


,,{g(o)£ 


e~]ks2dS  + 


(g(f)  ~  G( 0)) 

-2jks 


~ e~’ks2ds )  .  (3.22) 

ds  J 


Integrating  by  parts  once  more  (again  assuming  k  lossy),  yields 

I  _  ejkq(r,,<f>')  |  g(Q)  f°°  e-jks2ds+  (GM G(°.))e-jW 

l  Jia  2jksa 


jl  r d  [g(*)-c(o)i  ,-jk 

Jk  Jsa  ds 


^dsj  . 


(3.23) 


A  new  function  H(s)  is  defined  by 


and  written 


ff(.)  =  ±  [C(i)  -  C(0) 

ds  s 


H(s)  =  H (0)  +  (H(s)  -  H(0)) 


(3.24) 


(3.25) 


as  done  before  with  G(s).  Putting  this  into  the  integral  yields 

I  =  eJk<lir>J)  / G(0)  /°°e~jA:',2ds  +  GMjL^le-jk»a2 
l  Jsa  2  Jksa 

+ -u<  o)  r.-^ds + -l  r 

2jk  ’  JSa  2 jkJSa  [  -2 jks  J  ds  J 


(3.26) 


Integrating  by  parts  again  and  ignoring  the  unintegrated  term  yields 


Vj  ' -■  -VA ' *■  *  I-.*  f.-  /  /  vvA  -  Ov*  •  *.-  . * 


(3.34) 


V 

s' 


The  derivation  of  the  following  terms  is  shown  in  Appendix  A. 


yg 


-» 

tv 


H{sa)  = 


G(a0)  =  0 

saG'(sa)  +  G(0) 


Sa 4 


(3.35) 

(3.36) 


Using  these  relations  in  the  expression  yields  a  more  physically  revealing  form  of 
the  integrand  which  is  given  by 

H{  0) 


C(0)  + 


2jk 


ah) 


H{ 0)  (1  -  F(ksa2)) 


2  jk  2 jksa 


+e 


-ikr  gM 
{2jksa)2 


(3.37) 


Note  that  different  wave  components  can  be  identified  as  follows: 

^  —  ^re/  I tran 

where 


(3.38) 


ir„  =  ^ir.y)l/(_v4Sa)v/Z|G(0)  +  ^M 

(3.39) 

W  =  c-,k' 

G(0) 

2  jksa 

H(0){l-F(ksa2))' 
2jk  2 jksa 

(3-40) 

,  O'M 

d'<<  (2 jk..)* 

• 

(3.41) 

The  reflected  field  component,  Irf/,  containing  the  step  function,  corresponds 
to  the  surface  reflected  field.  It  vanishes  in  the  shadow  region  bounded  by  the 
RSB  and  exists  only  in  the  lit  region  beyond  the  RSB  when  sa  <  0.  The  higher 
order  term  is  required  because  the  first  order  term  vanishes  at  the  tip.  There  is  a 
tip-diffracted  or  end  point  term,  I diff-  This  term  is  the  same  term  that  was  found 


in  the  non-uniform  solution.  It  represents  the  field  diffracted  from  the  tip.  Lastly, 


there  is  the  transition  term,  Itran-  This  factor  keeps  the  total  solution  uniform 
across  the  RSB,  and  it  reduces  to  zero  in  the  deep  shadow  region.  Moreover,  the 
combination  of  all  three  components  provides  continuity  across  the  RSB. 

The  scattered  magnetic  field  is  then  given  by 


+  1 


tran 


+  ldi}}]  d<t> 


(3.42) 


or  more  specifically 


This  expression  will  not  yield  the  desired  expression  for  the  UTD  tip  diffracted 
field  as  yet.  One  problem  is  that,  in  this  form,  there  is  no  term  that  exists  on  the  lit 
side,  that  can  be  attributed  to  the  tip.  Before  the  modification  is  undertaken,  it  will 
be  shown  that  this  solution  is  uniform  across  the  reflection  boundary.  In  addition, 
the  behavior  far  from  the  reflection  shadow  boundary  in  the  deep  shadow  will  be 
investigated.  The  transition  term  plays  an  important  role  because,  in  general, 
the  reflection  boundary  transition  region  extends  into  almost  the  entire  region  in 
which  the  reflected  field  is  present.  Once  this  is  completed,  the  expression  will  be 
rewritten  and  the  UTD  tip  diffracted  field  identified.  It  will  be  shown  that  this 
new  form  yields  the  same  limits  as  the  quasi  UTD  form  near  the  reflection  shadow 
boundary,  and  far  from  the  reflection  shadow  boundary.  It  also  approximates  the 
quasi  UTD  solution  everywhere  within  the  transition  region. 

In  the  next  section,  the  expression  for  the  saddle  point  contribution  will  be 


determined.  It  will  be  shown  that  the  saddle  point  migrates  as  a  function  of  <f>\ 
and  the  implications  of  this  behavior  will  be  discussed. 

3.3  LOCATION  OF  THE  SADDLE  POINT 


The  saddle  point  rs  is  where  the  derivative  of  the  phase  function  in  the  inte¬ 
grand  of  (3.17)  goes  to  zero  where  the  phase  function  is  given  by 


q(r',<j>')  =  r'  cos  7  -  yr2  +  r2  -  2rr'  cosf 


(3.44) 


Then,  its  derivative  is 


d  ,  ,  r  —  r  cos  f 

—^{r  ,  <f> )  =  cos  7 - -=====  . 

r  y  r2  +  r2  —  2 rr'  cos  f 


(3.45) 


Setting  this  equal  to  zero  yields 


rs 


r(cos  f  ±  cos  7 


sinf 
sin  7 


(3.46) 


Substituting  this  result  into  Equation  (3.45),  it  can  be  seen  that  the  +  sign  satisfies 
the  original  equation;  therefore, 


/  sinf,  . 

rs  =  r(  cos  c  + cos  7- - ).  (3.47) 

sin  7 

Figures  4  thru  6  show  plots  of  the  migration  of  the  saddle  point  as  a  function  of 
<t>'  for  an  observer  located  in  specific  regions. 

In  the  lit  side  of  the  reflection  shadow  boundary,  the  saddle  point  rs  becomes 
positive  for  some  <f>'\  however,  it  is  always  negative  on  the  shadow  side  of  the 
reflection  shadow  boundary.  This  ties  in  to  the  existence  of  the  reflected  field. 


Recall  that 


igure  4:  Migration  of  Saddle  Point  with  <f>'  for  an  Observer  in  the  Deep  Shadow 

Side  of  the  Reflection  Boundary 


Figure  5:  Migration  of  Saddle  Point  with  <f>'  for  an  Observer  at  the  Reflection 
Shadow  Boundary  Corresponding  to  the  Specular  Line  of  Reflection 
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Figure  6:  Migration  of  Saddle  Point  with  <j>'  for  an  Observer  in  the  Deep  Lit  Side 

of  the  Reflection  Shadow  Boundary 

and  that  the  reflected  field  has  a  unit  step  function,  U{  —  \/ksa).  When  rs  <  0, 
then  sa  >  0  and  there  is  no  reflected  field.  This  follows  since  rs  <  0  occurs  in  the 
shadow  side  of  the  reflection  shadow  boundary.  When  r9  >  0,  then  sa  <  0,  and 
there  can  be  a  reflected  field  component.  It  is  worth  noting  that  even  in  the  lit  side 
of  the  reflection  shadow  boundary,  rs  is  not  always  greater  than  zero,  and  that  it 
varies  as  a  function  of  <f>'.  That  is,  the  reflected  field  exists  only  for  certain  aspects. 
Also,  the  saddle  point  approaches  the  end  point  which  therefore  requires  a  more 
careful  treatment  of  the  integral  in  (3.17).  That  is  the  reason  why  the  asymptotic 
approximation  to  the  integration  in  (3.17)  was  performed  as  shown  in  the  previous 
section. 

This  saddle  point  migration  was  also  observed  by  Molinet  [39]  in  his  research. 

In  the  next  section,  the  behavior  of  the  integral  at  the  reflection  shadow  bound¬ 
ary  and  in  the  deep  shadow  side  of  the  reflection  boundary  will  be  investigated. 
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3.4  LIMITING  CASES 


*W»!« 


Two  limiting  cases  are  investigated  for  the  result  in  (3.37).  Near  the  RSB, 
and  in  the  deep  shadow  side  of  the  reflection  shadow  boundary. 

Substituting  the  actual  expression  for  H  from  (3.24)  into  the  integral  in  (3.34), 
and  recalling  that  G(s0)  =  0,  yields  the  following  in  <f>': 

I M\fc)  ^  G(0)  +  ^^ 

+e-i*r  m  _L_\1  _  G"(°)  (1  ~/(fcjg!))1 

[2jksa  j  V1  2jksa2)\  4 jk  2jksa 

<3<9> 

This  is  the  expression  from  which  uniformity  is  proven.  The  small  and  large 
argument  forms  of  the  transition  function  F(X)  are  given  in  [19]  as, 


F[X) 


1  -  Jjx  for  X>>1,  and 
yJnjXei*  for  X<<1. 


(3.50) 


Near  the  shadow  boundary,  X  — ►  0;  therefore,  the  small  argument  form  of 
the  transition  function  is  used.  In  terms  of  specific  components,  the  expression 


becomes 


lr(f  =  e’W^Ui-Vksa)^  G(0)  +  ^J^  (3.51) 

G»(0)  [l  ~ 


2  jksa 


,  _  CM 

d'"  (2 jkstf 


(3.52) 


(3.53) 


V-VAV..'.','.  A.  A.  .%  A  A.  A  A.  A  A-  A-  A  A. 
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Next,  G  is  expanded  in  a  Taylor  series  about  sa.  This  is  given  by 


C(>)  =  G(sa)  +  (3-  sJG'M  +  (*  ^-a"M 

(3.54) 

from  which 

G  (s)  —  G  (sa)  +  (s  sa)G  (sa)  +  ^  G  (sa)  +  •  •  • 

(3.55) 

G"(s)  =  G"(sa)  +  (s-  Sg)G'"(Sg)  +  . . . 

(3.56) 

Recall  that  G(sa)  =  0;  therefore,  evaluating  at  s  =  0  yields 

.  2 

G(0)  =  ~SgG' (sa)  -f-  ——Gri(sg)  +  .  .  . 

(3.57) 

G'(0)  =  G'{Sa)-SgG"(Sg)+... 

(3.58) 

G"(0)  =  G"(Sa)  +  ... 

(3.59) 

Substituting  these  into  Equations  (3.51)  and  noting  that  sa2  -  r  =  9(0)  -f  sa2  = 
q(rf),  yields 


rG"(sa) 
4  jk 

+  e? 


'rtf 

Itran 

idtff 


=  e 


+  o(«a) 


=  e~’k 


^  1  ■ 


rG'(sa) 


L  2 jk  (2j'fcsa)2  J 


(3.60) 


—  ikr  G'(sa) 

(2 jksa)2' 


(3.61) 

(3.62) 


Approaching  from  the  shadow  side  of  the  reflection  shadow  boundary,  sa  — *  0+ 
(sa  >  0),  the  integrand  becomes 


I(r,  <t>,  k)  =  hran  +  Idiff 

=  iejkqM)  [*~  G"{Sg)  lkrG'{Sg 

2  \  jk  4jk  2 jk 


(3.63) 

(3.64) 


Approaching  from  the  lit  side  of  the  reflection  shadow  boundary,  sa  — ►  0  (sa  <  0), 
the  integrand  becomes 

!(»■,<£,£)  =  Iref +  hran +ldxff 

=  cjkqir.,4?)  G"(sa)  /  l\  -jhr&M 
\jk  4 jk  V  2/  '  2jk 

2  4jA  2jk  1  1 

This  requires  that  the  limit  as  «a  -*  0  be  taken,  which  means  that  — »  0. 
This  occurs  at  the  reflection  shadow  boundary,  where  f  +  7  =  it.  In  addition, 
9(rfi</>/)  =  —  r  at  the  reflection  shadow  boundary 

From  the  above,  it  has  been  shown  that  the  integrand  in  <f>\  and  hence  the  in¬ 
tegral  over  <t>\  is  continuous  across  the  reflection  shadow  boundary  and  the  solution 
is  uniform. 

Far  from  the  RSB  the  large  argument  form  of  the  transition  function  is  used. 
The  behavior  far  from  the  RSB  on  the  shadow  side  of  the  RSB  will  be  investigated. 
In  this  case,  one  finds  that 


I  ref  =  0 

(3.66) 

liran  =  0 

(3.67) 

T  —ikr  G'{Sa) 

ld'<>  '  (2 jk..r 

(3.68) 

The  above  results  follow  because  on  the  shadow  side  of  the  reflection  bound¬ 
ary,  sa  >  0  and  the  step  function  term  containing  the  reflected  field  vanishes 
because  there  is  no  reflected  field  in  this  region.  In  the  transition  term,  there  are 
1  -  F(ksa 2)  factors.  Using  the  large  argument  form  of  the  transition  function, 
these  asymptotically  go  to  zero.  Therefore,  I  reduces  to  the  tip-diffracted  term 
alone  such  that 


I(r,<M)  =  I  dtff- 


(3.69) 


As  previously  mentioned,  once  integrated,  this  yields  the  same  value  for  tip  diffracted 
fields  as  the  exact  evaluation  of  the  far  zone  PO  integral.  Note  that  on  the  lit  side 
of  the  RSB,  the  argument  of  the  transition  function  F  is  never  sufficiently  large; 
thus,  unlike  the  shadow  side,  the  transition  region  extends  over  the  entire  lit  side 
of  the  RSB. 

3.5  UTD  TIP  DIFFRACTED  FIELD 

The  UTD  tip  diffracted  field  contains  a  factor  (transition  function)  which 
keeps  it  bounded  at  the  RSB.  In  the  quasi  UTD  form  from  the  previous  section,  the 
singularities  at  the  GO  boundaries  were  corrected  additively.  This  is  similar  to  the 
way  that  the  relevant  singularities  cancel  at  the  RSB  in  the  UAT  form  described  in 
the  beginning  of  the  chapter.  Molinet  [39],  has  independently  worked  the  problem 
of  scattering  from  a  cone  using  the  PO  approximation.  This  was  done  for  the 
special  case  of  axial  incidence,  and  he  basically  writes  his  high  frequency  solution 
in  the  form  of  the  uniform  asymptotic  theory  (UAT)  [40]  which  follows  from  the 
use  of  asymptotic  methods  discussed  in  Felsen  and  Marcuvitz  [36].  The  additive 
correction  in  the  UAT  form  requires  special  treatment  to  overcome  the  numerical 
difficulties  encountered  in  programming  this  type  of  solution.  The  current  UTD 
expression  does  not  encounter  this  difficulty.  In  addition,  it  enables  the  isolation 
of  a  uniform  tip  diffracted  field.  The  tip  diffracted  field  is  extremely  useful  for 
both  the  scattering  and  the  antenna  problems. 

The  UTD  form  for  the  scattered  field  is  given  by 

U,c  «  V9r°  +  (3.70) 

where,  for  the  cone,  the  first  term  is  the  GO  reflected  field  component,  and  the 
second  term  is  the  tip  diffracted  field  component.  It  can  be  seen  that  in  order  to 


get  the  expressions  from  the  previous  section  in  this  form,  one  must  first  recognize 


from  (3.70)  and  (3.39)  that 


u  r  =  i. 


=  G( 0) 


m 


(3.71) 


-  W  +  U//  =  .-'*'[^[r(W)-(.-^)]  (3.72) 


fl(0)(l-f(*»»2))j  (  r,kr  gM 
2  jk  2jksa  (2 jksa)2 


(3.73) 


where  it  was  previously  given  that 


=  .-I*  f  cm  rrrt.  2,  /,  i  \\  _  m  (i  -  n^)) 

2jksa  a  \  2jfcsa2/J  2 jk  2jksa 


-  „-jkr  ^(-Sa) 


(3.74) 

(3.75) 


ld'H  “  (2^'  (3  75) 

One  can  futher  recognize  that  the  terms  in  Itran  are  important  only  near  the 
reflection  shadow  boundary;  whereas,  they  are  negligible  away  from  the  reflection 
boundary.  Near  the  boundary,  sa  is  small  and  the  Taylor  series  expansion  for  G(0) 
and  H(0)  can  be  used. 

Having  the  expressions  in  terms  of  the  end  point,  rather  than  the  saddle  point 
also  gives  an  even  clearer  picture  of  the  mechanisms  involved.  The  scattered  field 
will  have  two  identifiable  sources;  the  end  point  (tip)  diffracted  field,  associated 
with  the  UTD  tip  diffracted  field,  and  the  surface  reflected  field.  Once  the  expres¬ 
sions  are  derived,  it  will  be  shown  that  they  reduce  to  the  same  limits  as  the  quasi 
UTD  expression  from  the  previous  section. 


Recall  that  H( 0)  =  ^G"(0)  and  that 


=  2 

G(0)  =  —  saG'(sa)  H — — Gw(sa)  +  . . . 
G'(0)  =  G'{sa)  -  saG"{sa)  +  ... 

G"{  0)  =  Gn(sa)  +  — 


(3.76) 

(3.77) 

(3.78) 
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Substituting  this  into  the  expression  for  Itran  in  (3.74)  yields 


Ifran  '  e 


—jkr  \ G^[Sa)F{kSa 


[  4jk  2 jks 

-■jkr  G'(sa) 

(2 jksa)2 

The  last  term  cancels  therefore, 

Ud  =  Itran  +  I*//  =  e">fcr 


(3.79) 


G"(«a)  F{ksa2) 


+ 


(3.80) 


[  4 jk  2jksa 

Once  again,  the  small  and  large  argument  forms  of  the  transition  function  F(X) 
are  given  below: 

f  1  -  2 jx  f°T  X>>1,  and 


F(X)  « 

{  y/njXe^  for  X<<1. 
Far  from  the  reflection  shadow  boundary, 


(3.81) 


Hk’S)  *  '  -  (3  82) 

Assuming  the  observation  point  is  on  the  shadow  side,  then  the  reflected  field, 
U 9r° ,  is  zero,  and  all  that  remains  is  \Jd.  Inserting  the  large  argument  form  shown 
above  into  yields 


U“  = 


p-jkr 

1 

1 

1 

1 

\  (CM 

C 

4jk 

.2 jksa  {2jk)2sa3\ 

T 

2 jksa2 

\  V  2 jk 

s 

4Jk  ). 

t~]kr 

'  G'{sa) 
(2  jksa)12 

i  +  o{k  3) 

«  e 


-jkr  G'(Sg) 
(2 jksa)2 


(3.83) 


which  is  the  same  limit  as  in  the  previous  section.  As  before,  this  is  also  the  value 
for  the  tip  diffracted  fields  that  the  exact  evaluation  of  far  zone  PO  integral  yields 
in  this  region. 


Close  to  the  reflection  shadow  boundary,  the  small  argument  of  the  transition 
function  is  used.  That  is 


F(ksa2)  ss  \/njk\sa\e}kSai 


(3.84) 


Hence,  for  ka— ♦  0 


1  -  F(ksa2)  «  1 


Therefore,  noting  that  sa2  -  r  =  9(0)  +  sa2  =  q{rs),  then 


,  -ikrG>{*a) 


(3.85) 

(3.86) 


(3.87) 

(3.88) 


Ur  «  e^\[^^-U(-Vksa) 

V  JK  “tj  K 


(3.89) 


Approaching  from  the  shadow  side  of  the  reflection  shadow  boundary,  sa  — ►  0+ 
(sa  >  0),  the  integrand  becomes 


Usc  =  Uc 


_  ijkqM)  [*’&{**)  1  c-ikr&M 
2  V  Jk  4jk  +  2  jk 


(3.90) 

(3.91) 


Approaching  from  the  lit  side  of  the  reflection  shadow  boundary,  sa  -+  0  (sa  <  0), 
the  integrand  becomes 

usc  =  ur + vd 

=  gjkqM')  r*~G"(*a)  (l  ]\  ■  -jfcrgM 
y  jk  4  jk  V  V  2jk 

=  [Z<lM  +  e-ytrgM  (3.92) 

2  V  4jfc  2jfc 


Again,  this  is  the  same  limit  as  before.  Also,  as  before,  this  requires  that  the  limit 
as  sa  — ►  0  be  taken,  which  means  that  r3  — *  0.  This  occurs  at  the  reflection  shadow 
boundary,  where  f  +  qr  =  n.  In  addition,  q[rs,<f>')  =  —r  at  the  reflection  shadow 
boundary 

Investigation  of  shows  that  it  is  a  function  of  only  the  end  point.  For  unit 
amplitude  plane  wave  incidence,  using  the  tip  as  the  phase  reference,  l/*(tip)  =  l. 
Therefore,  letting  the  point  Q  be  the  tip,  yields  the  usual  form  given  by 

\]d  =  U*(Q)  DTe~}kr  (3.93) 


where 

— <**> 
One  should  recall  that  G'  and  G"  are  vector  quantities,  but  for  convenience,  were 
not  written  in  vector  notation;  therefore,  is  a  vector  quantity.  This  is  used  to 
determine  the  contribution  from  tip-base  interaction  for  the  finite  cone.  It  can  also 
be  used,  via  reciprocity,  to  determine  the  radiation  from  an  infinitesimal  current 
element  on  the  cone  surface  due  to  this  interaction. 


The  scattered  magnetic  field  is  then  given  by 


hs  ~  J2*  [ur + ud]  d<f>' 


Moreover,  The  magnetic  field  from  the  tip  is  then  given  by 


Hrip *  J2*  DTt~]kr d<f> 


(3.95) 


(3.96) 


where  Dj  was  given  in  (3.94). 


3.6  RESULTS 


In  this  section,  plots  of  the  scattered  field  for  various  cone  configurations  will 
be  shown.  In  addition,  the  uniform  PO  solution  is  compared  to  the  far  zone  PO 
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solution  for  the  scattered  field  from  Chapter  II.  Both  will  be  compared  to  Felsen’s 
small  angle  cone  result  [5]. 

Theory  is  compared  to  experiment.  The  tip-base  interaction  shown  in  Figure 
19  is  included  in  the  computation  of  the  RCS  for  a  finite  flat-backed  in  the  result 
given  by  Burnside  and  Peters  [10],  and  it  is  shown  that  this  mechanism  corrects 
the  previous  discrepancy  between  calculation  and  measurement. 

Let  us  begin  with  the  comparisons.  The  H £  component  found  in  terms  of 
both  the  uniform  and  exact  far  zone  PO  solutions  are  plotted  in  Figures  7  thru 
18  with  Felsen’s  small  angle  cone  result.  There  is  reasonable  agreement  over  most 
of  the  shadow  side  of  the  reflection  boundary  between  Felsen’s  result  [5]  and  the 
uniform  PO  solution  for  the  case  of  axial  incidence.  The  agreement  is  excellent 
for  small  cone  angles  where  Felsen’s  result  [5]  is  the  most  accurate;  however,  the 
agreement  holds  up  even  for  cone  angles  up  to  15°.  When  the  incidence  is  nonaxial, 
the  agreement  between  the  present  result  and  the  one  obtained  using  the  result  in 
[5]  is  still  good  for  the  narrow  angle  cone,  but  is  1  or  2  dB  worse  than  that  for 
axial  incidence  in  the  case  of  wider  angle  cones.  The  present  uniform  PO  solution 
and  the  far  zone  PO  solution  are  valid  when  the  cone  is  fully  illuminated.  This  is 
true  when  the  incidence  angle  is  smaller  than  the  internal  cone  angle.  Numerical 
results  were  also  obtained  for  incidence  angles  that  were  at  or  near  this  limit  to 
see  how  the  uniform  solution  behaved.  As  long  as  one  is  a  few  degrees  away  from 
this  boundary,  the  results  look  fine  compared  to  those  far  from  this  boundary; 
however,  when  evaluated  right  on  this  boundary,  as  shown  in  Figure  12,  the  curve 
shows  an  unexpected  drop  of  approximately  8  dB  at  the  cone  surface.  One  would 
not  expect  the  field  to  vary  this  much  in  the  span  of  1°;  therefore,  it  is  assumed 
that  this  is  a  numerical  problem;  however,  this  could  also  be  a  result  of  grazing 
incidence  when  the  incident  field  lies  along  this  boundary. 
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Figure  9:  Comparison  of  Hi  for  a  =  15°,  0  =  0°,  <f>  —  O' 
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Figure  13:  Comparison  of  Hi  for  q  =  15°,  0  =  0°,  <j>  —  45c 


The  UTD  tip  diffracted  field  is  used  to  determine  the  tip-base  interaction 
shown  in  Figure  19.  The  tip-base  interaction  is  included  in  the  GTD  analysis 
conducted  in  [lOj.  In  Figure  20,  the  plots  of  Burnside  and  Peter’s  [10]  GTD 
solution  for  the  RCS  of  the  flat-backed  cone  are  compared  to  measurement.  In 
Figure  21,  the  GTD  solution  in  [10]  is  compared  with  a  moment  method  (MM) 
solution  for  the  finite  flat-backed  right  circular  cone.  It  is  noted  that  the  MM 
solution  shows  essentially  the  same  behavior  as  the  measurements;  furthermore, 
the  GTD  solution  in  [10]  is  in  error,  with  respect  to  both  measurements  and  the 
MM  solution,  by  approximately  3  dB  over  the  entire  frequency  range  8  <  ka  <  12. 
In  Figure  22,  the  comparison  of  the  GTD  solution  in  [10]  which  is  modified  to 
include  the  tip-base  interaction  based  on  the  current  work  is  shown;  indeed,  it  is 
evident  from  this  figure  that  the  inclusion  of  the  tip-base  interaction  essentially 
removes  the  error  which  results  if  the  latter  interaction  is  not  included.  The  ripple 
in  the  MM  solution  is  probably  caused  by  the  interference  between  the  direct  tip 
and  base  diffraction  contributions.  It  is  noted  that  the  direct  diffraction  by  the  tip 
is  not  included  in  Figures  20  and  21  or  in  the  modification  of  the  GTD  result  of 
[10]  shown  in  Figure  22. 


CHAPTER  IV 


VECTOR  WAVE  FUNCTION  SOLUTION 

4.1  INTRODUCTION 

In  this  chapter,  the  eigenfunction  expansion  of  the  dyadic  Green’s  function 
for  the  semi-infinite  cone,  derived  in  Appendix  C,  is  used  to  describe  the  fields  of 
an  arbitrary  electric  point  source.  The  expansion  is  then  specialized  to  the  case 
where  an  electric  current  point  source  is  located  on  the  forward  axis  of  the  cone. 
The  resulting  expressions  are  investigated  as  the  observation  point  nears  the  tip. 
Lastly,  the  problem  of  the  radiation  by  circumferential  slot  antennas  on  the  cone 
is  solved. 

4.2  SOLUTION  USING  THE  DYADIC  GREEN’S  FUNCTION 

The  geometry  for  the  semi-infinite  cone  is  shown  in  Figure  23. 

Beginning  with  expression  (C.34)  in  Appendix  C,  where  E  and  H  are  the  total 
electric  and  magnetic  fields,  and  T  is  the  dyadic  Green’s  function,  and  substituting 
the  expression  for  the  equivalent  electric  surface  current,  J8,  which  is  given  in 
terms  of  the  magnetic  field  at  the  surface  by  J^R')  =  n1  x  H(R')  yields 

R  iToUn  V  T  }  =  fv  -MR'>  •  m'W 
+juv  j  n'  x  H(R')  •  F(R',R)dS'  +  f  K«(R')  •  (V'  x  f  (R',R ))dS'  (4.1) 

where  V  is  the  volume  enclosed  by  the  surface  S,  and  the  vector  R'  is  the  source 
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Figure  23:  Geometry  for  Semi-Infinite  Cone 
vector  given  by  R'  =  [R' ,0' ,<t>')  and  R  is  the  observation  vector  given  by  R  = 

(RJ,4>). 

The  Green’s  function  has  been  chosen  such  that 

x  P(R',»)]  =  -*  x  F(R',R)]  =  0.  (4.2) 

The  eigenfunction  expansion  of  the  dyadic  Green’s  function  for  the  cone,  in  terms 
of  the  vector  wave  functions,  which  meets  these  conditions  is  derived  in  Appendix 
C.  It  is  given  by 

'  ^  nJ,)(R)nJ4)(R')  „  „ 

f(R  R'l  =  - nT - +  ^  - TT -  R<R  U3) 

’  v  mJ4)(R)mJ,)(R')  r  nJ4)(R)n11}(R')  ^  , 

Lp  -rip - +  ^9  — rf -  R>  R 

Since  p  and  q  are  distinct,  application  of  Equation  (4.2)  yields 

e  x  M^R)]^  =  i>  x  Np»(R)]  =  0.  (4.4) 


where  p  =  %  mu  and  q  =  %  mi-  Evaluating  these,  it  is  seen  that 


**mp1|(R)]<o  =  o 

(4.5) 

which  requires 

4fZ"(  co.»)|  =0, 

J  6=6q 

(4.6) 

and 

*xN<V)]^o=° 

(4.7) 

which  requires 

^^(cos  0q )  -  0. 

(4.8) 

Commuting  the  dot  and  cross  products  in  the  first  surface  integral  yields 


n'  x  H(R')  •  r(R', R)  =  -H(R')  •  n'  x  T(R',R),  (4.9) 

but  from  Equation  (4.2),  n'  x  r(R\R)  =  0;  therefore,  the  first  surface  integral 
goes  to  zero.  In  other  words,  the  Green’s  function  accounts  for  the  presence  of 
the  cone.  Moreover,  the  cone  is  assumed  to  be  a  perfect  conductor;  therefore,  the 
equivalent  magnetic  surface  current,  Ks,  is  given  by  Ke  =  0.  This  eliminates  the 
second  surface  integral. 

Once  these  preliminary  steps  are  carried  out,  all  that  remains  is 

E(R)  =  jkZc  f  JV(R')  •  f(R',R)dK'.  (4.10) 

JV 

This  can  be  rewritten  in  a  different  form.  It  can  be  shown  that  for  a  vector-dyadic 
pair 

a -7  =  7-a 

Letting  a  =  JvfR')  and  /  =  r(R',R),  and  recalling  that 


(4.11) 


the  electric  field  can  then  be  written  as 


E(R)  =  jkZc  f  r(R,R')  •  Jv(R')^' 

JV 


(4.13) 


where  w/x  --  kZc ,  k  =  Uy/jie,  and  Zc  —  \J~^' 

Using  R"  as  a  dummy  variable  of  integration,  and  defining  JV(RW)  in  terms 
of  a  point  electric  current  source  at  R',  the  volume  current  source  can  be  written 


as 


JV(R")  =  pe(R")*(R"  -  Re¬ 

inserting  this  expression  into  Equation  (4.13)  with  Rw  yields 

E(R)  =  jkZc  f  f(R,R")  pe(R"WR"  -  R -)dV" 
Jv 

Applying  the  sifting  property  of  the  delta  function  yields 


(4.14) 


(4.15) 


E(R)  =j*Zcr(R,R')-pe(R') 


(4.16) 


which  yields  an  expression  for  the  electric  field  in  terms  of  a  point  current  source 
of  strength  pe  at  R'. 

Once  the  E  field  has  been  determined,  it  is  an  easy  task  to  determine  the  H 
field  and  then  the  surface  currents  that  will  be  present  on  the  cone.  For  the  sake 
of  completeness,  it  is  noted  that 

H(R)  =  -V  x  f  (R,R')  •  Pe(R').  (4.17) 

4.3  EXPLICIT  EXPRESSIONS  FOR  THE  FIELDS  OF  AN  ELEC¬ 
TRIC  CURRENT  POINT  SOURCE  IN  THE  PRESENCE  OF 
THE  SEMI-INFINITE  CONE 


It  is  still  assumed  that  the  plane  of  incidence  is  the  x-z  plane.  This  can  be 
done  without  any  loss  of  generality  since  the  cone  geometry  is  very  symmetric. 


Figure  24:  Geometry  for  Application  of  Dipole  Moment 


If  another  aspect  were  desired,  one  could  just  rotate  the  coordinate  system.  The 
choice  of  the  x-z  plane  as  the  plane  of  incidence  simplifies  the  expressions. 

In  general,  the  strength  of  the  electric  current  point  source  can  be  written  as 


Pe  =  Pet' 


(4.18) 


Two  orientations  of  the  current  will  be  dealt  with.  They  are  e'  =  O'  and  e1  =  <(>'. 
Any  arbitrary  polarization  of  the  incident  field  arriving  from  a  distant  source  can 
be  expressed  as  a  linear  combination  of  these.  Figure  24  shows  the  geometry  and 
the  location  of  the  point  source. 

The  point  source  is  placed  such  that  R'  =  d  >  R\  therefore,  the  representation 
of  the  Green’s  function  that  is  valid  for  R  <  R'  is  used  because,  eventually,  the  field 
near  the  cone  tip  is  desired.  The  source  position  vector  is  given  by  R/  =  (d.,/3,0) 
since  the  plane  of  incidence  is  assumed  to  be  the  x-z  plane.  The  incidence  angle 
0  is  chosen  such  that  O<0<Oo  =  n-a,  where  a  is  the  cone  half-angle.  The 
selection  of  the  required  Oq  or  a  is  needed  to  solve  the  eigenvalue  problem  for  u 


% 


The  fields  due  to  point  sources  which  are  oriented  transverse  to  R!  can  be 
found  by  letting  pe  =  pe6'  and  pe  =  pe4>'  and  by  letting  pe  =  1  for  convenience. 
The  6'  and  <f>'  in  the  superscripts  show  which  point  source  orientation  generates 
the  corresponding  fields  and  currents.  This  yields 


Ee,(R)  =  jkZcT(R,R')  •  O' 

(4.19) 

and 

E*'(R)  =jkZcT{RyR') 

(4.20) 

Each  case  separately  is  worked  separately.  In  the  O'  case 

E*'(R)  =  jkZc 

'—mJVjmJV')  ^  y'N51|(R)N5,|(R') 

■  P  9 

(4,0,0) 

(4.21) 

This  yields 

mJ,4)(r')  • 

Al  _ PU*{c°S  /?)  (2)  t x„(4)  z, 

<W,0)-m  sin/?  K  td  ~  M°™'(R  )  ' # 

(4.22) 

and 

N^4)(R')  -S' 

(</,/>, 0)  ki  <10  ^  lcosW )  0 

(4.23) 

so  for  the  6'  polarization 

p  =  omu 

(4.24) 

q  -  em 7 

(4.25) 

Therefore,  the  electric  field  in  this  case  is  given  by 


E*(R)  =  jkZc 


£  2)(W)mP™(cM^ 


L  P 


n£ 


sin/? 


E*'(R)  is  found  in  a  similar  manner.  Applying  Equation  (4.20)  yields 


E*  (R)  =  jkZc 


E 

p 


Mp  (R)Mp  (R') 


nt 


*'+E 

9 


Ng  (R)N^  (R;) 


4>' 


Performing  the  indicated  dot  products  yields 


WO,0)  =  ■k*l(“)^p™(cos/?)  =  mS41(R') -i> 


and 


N’  (R  > '  *  1  M«i "  ~SnF - id - -  N-"->(R  > '  * 

so  for  the  4>'  polarization 

p  =  emu 
q  =  om-y 

Therefore,  the  electric  field  in  this  case  is  given  by 

E*'(R)  =  jkZc 


(d,p,  0) 

(4.27) 


(4.28) 

(4.29) 

(4.30) 

(4.31) 


+E 

9 


Ni^(R)  (kdh\i](kd))'  P?{ cos/?) 

-  — m — - - 


.(2)/ 


kd 


sin/? 


(4.32) 


Therefore,  it  is  seen  that  the  location  and  direction  of  the  sources  determines 
the  modes  which  are  summed.  As  seen  in  the  E^  case,  the  non-zero  values  of 
M  •  0'  and  N  •  O'  determined  the  specific  behavior  of  p  and  q.  As  was  said  earlier, 
the  choice  of  the  plane  of  incidence  made  the  expressions  simpler;  however,  no 
generality  was  lost  as  a  result. 

Equations  (4.26)  and  (4.32)  contain  similar  expressions  involving  the  Legendre 
functions.  To  simplify,  write 

P?(  cos*) 


S(a,m,0)  =  m 


sin  0 


(4.33) 
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and 


T(a,m,t)  =  -P”(cos»). 


(4.34) 


These  can  be  inserted  into  the  previous  equations  for  E^  (R)  and  E^  (R). 
Using  recursion  formulas  found  in  [18],  S  and  T  can  be  written  as 

p-+1(c°s  0)  +  (a  -  m  +  l)(ct  +  m)P™~1  (cos  0) 


S(a,m,8)  = 


2  cos  8 


(4.35) 


and 


T(a,m,0)  = 


P™  1  (cos  0)  (a  -  m  +  1)(q  +  m)  -  P™+1(cosO) 


(4.36) 


Recalling  the  relationship  between  M  and  N,  and  applying  Maxwell's  equa¬ 
tions  to  the  expressions  for  the  E  fields  yields 


H*'(R)  =  -k 


EN[mi/(R),(2)a  ,,c, 

- - — hi  (kd)S(i/,m,0) 

L  P  v 


9 


mSH,(R)  (kdh^(kd))' 


(2), 


kd 


T(~i,m,0) 


(4.37) 


H*'(R)  =  -k 


Y,  N%(R-)-(-/.!,2>M)rKm.g) 


+ 


E 


M{X(R)  (kdh\£,(kd)Y 


(2)  i 


kd 


S(~i,m,0) 


(4.38) 


Evaluating  the  H  fields  on  the  surface  of  the  cone  at  Rq  =  (Ro,0o,4>),  and 
calculating  the  surface  current  J8(Rq)  given  by  n  x  H(Rq)  where  n  =  —0,  yields 


jJ(Ro)  =  * 


9  X  Nomi/(Ro)  l(2)/,  j\  cl  a\ 
- - hi  '(fcd)S(i/,m,/?) 

Up 


i  p 


,  ^8  xm[]1,(R0)  (kdh^(kd))'^ 

+  E - t: - kd - 


(4.39) 
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Jj(»o)  =  k 


E0  x  Nemi/(Ro) ,  L(2)n  jut/ 

- - - (-hi '( kd))T(i/,m,P ) 


L  P 


n. 


t  y  >  *  Ml^(Ro)  (UhSPMY  m  „ 


kd 


(4.40) 


where 


0  X  nSL(Ro) 

- d>i/(u  +  l ) P? (cos *6) sin m4> 
kRo 

(4.41) 

9  x  M^(R0) 

=  -R^zP?l(cos9o)j~l(kRQ)cosm<j> 
du  ' 

(4.42) 

9  x  nSIL(Ro) 

-  f  _  .  Ji/(kRo)  1  P™(cos<?o)  .  , 

-  kRo  }m  sin,0 

-  iffl'  +  1)  [  n  ^  P™(cos  0q)  cos  m<p 

kRo 

(4.43) 

9  x  M^(R0) 

=  -R~P™(c°sOo)h(kRo)sinm<j). 

(4.44) 

From  Appendix  B, 

_  27 r 

flmj/  =  ~T  "F 
JKt-m 

Tmi  =  jktrrC1^1  + 

where  the  normalization  integral,  W,  is  given  by 


(4.45) 

(4.46) 


cos  0)]2  sin0d0. 


(4.47) 


Substituting  Equations  (4.41)-(4.46)  into  Equations  (4.39)  and  (4.40)  the 


component  surface  currents  are  given  by 


jfa(*  o)  = 


jfe2  y- 

2tt  ^  i/(t/  +  l)Jm„ 


f  S(i/,m,/?)  cosm<£  1  (”  //.«  v  Jv(kRo)]  mP™(cos90) 

\  }  \MkR0>  +  sin*„ 


Etmj-i{kRo)  (  T(i/,m,/3)  cos  m<£  1  /  ,  (2)# 
'y(7  +  l)/m7\  S(Tf,m,/?)sinm0  J  |  *> 


(m)  +  ^lm^1sp”(cos#o) 


(4.4«) 


y'(Ro)  =  3—  E  -r2-  {  _T hmx)smmt  )  °*«o) 

u;  2*  ^  Jm„  \  T(i/,m,/?)cosm<£  J  "  v  ’  ^Rq  v  v  u' 


(4.49) 


The  expressions  for  the  surface  current  can  vary  as  either  a  function  of  the  az¬ 


imuthal  angle  <f>,  or  the  radial  distance  from  the  cone  tip  Rq. 


In  the  next  section,  the  case  for  which  the  dipole  moments  are  on  the  cone 


axis  are  shown;  in  other  words,  for  social  incidence.  The  recursion  formulas  will  be 


especially  useful  in  this  case. 


4.4  SPECIALIZATION  TO  THE  CASE  OF  A  POINT  SOURCE  ON 
THE  CONE  AXIS 


In  this  section,  the  surface  currents  for  the  case  when  the  electric  current 


point  source  is  on  the  cone  axis  are  determined.  The  source  on  the  cone  axis  is 


chosen  to  be  oriented  transverse  to  the  cone  axis  so  it  could  have  9'  and  </>',  or  x 


and  y  components.  In  this  case  0  —  0.  We  will  be  starting  with  Equations  (4.48) 
and  (4.49)  for  the  currents.  The  value  of  m  can  be  determined  by  investigating 
Equations  (4.35)  and  (4.36)  for  0  =  0.  In  both  cases  note  that 


PM  =  l 


(4.50) 


mm 


-  '  v  -  -  - 


this  implies  that 


P™(1)  =  0  m  >  0. 


(4.51) 


Looking  at  Equations  (4.35)  and  (4.36),  both  S  and  T  have  a  term  of  the  form 
P™_1(1).  Cleary  this  will  be  zero  when  m  -  1  >  0,  or  m  >  1.  Therefore, 

S(a,m,0)  =  T(a,m, 0)  =  j  2a(aQ+  *)  ™  =  J  (4.52) 

The  following  representations  for  the  current  due  to  an  axially  positioned  point 
source  are  given  by 


J^(R0)  =  k 


flu,  2 


Ox  m[^(R0)  (Jcdh™<JcW  +  l) 
V  rl7  kd  2 


.(2), 


(4.53) 


and 


J.Vo)  =  k 


E 


**  lv  2 


v-  9  x  M'li.fRo)  (kdh^’(kd))'  + 

+  V  Tn  M  2 

Evaluating  the  previous  expressions  for  m  =  1  yields 


.(2)/ 


(4.54) 


ftiz,  =  +  l  )hu 

jk 

Tn  =  + 

from  which 

+ 1)  1  _  jk 

2  Hi,,  2nliu 

Th  + 1)  i  _  i* 

2  Th"  2*7,, 
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(4.55) 

(4.56) 


(4.57) 

(4.58) 


where  /ja  is  given  as  in  Appendix  B  as 


rv  o 

ha  =  /  [Pa(cos0)]2sin0<f0.  (4.59) 

Jo 

For  v, 

Jq  [Puicos6)\2s'in0de  =  ^~~Pu{^se0)-^^P^{co8e0)  (4.60) 

and  for  7, 

J  \P'(cos6Q)}2sin0d0  =  ^  P^j  (cos  0O)~P^  {cos  0O).  (4.61) 

Recall  from  Appendix  D,  that  for  m=l 

d2 

U(*o)  =  J^-0P^os0o)  (4-62) 

L\(0O)  =  ~~P]{  cos<?0)  (4.63) 


When  these  functions  are  evaluated  on  the  cone  surface,  there  will  be  some 


factors  that  cancel.  Prior  to  the  cancellation,  the  component  expressions  for  the 
current  are  given  by 


where 


Vr  = 


9 

W  (p\  _  jJL  /  cos  <j>  1 
JsR\H-o)  “  2jr  \  sin*  J 

(4.64) 

0f 

t t>  \  _  /  -  sin  <t> 

-  2?r  \  cost 

(4.65) 

{***> +’?*»] 

1  pl  (cos  0O) 
sin  9q 

,(2)/ 


Y,  ~  |h!,2)  (fed)  +  —  ~~~  |  h(kR0)~P^(cos  0O) 


(4.66) 


and 


**  =  £ 


After  cancellation, 


)  u(2)n.j\^kR^  dU  a\ 

hi/  {kd)  Pv{™s$o). 

(4.67) 

~J^  (  COS  <j>  1  1 

2tt  sin  0O  \  sin  4>  j  R 

(4.68) 

)-  ~jk2  A 

0  27r  sin  0O  ^ 

(4.69) 

where 


V  <2"  +  >)  '>!•  (*'')  I  .  „  ,  .  MkR o)  1 
A*=  £-*»*„  «<«,) 

1  lo* 


(4.70) 

h[?]  {kd)  j^kRo) 

Ll(0o)  kRo  • 

(4.71) 

„  =  (-» 
The  source  can  be  placed  in  the  far  field.  This  allows  the  use  of  the  large 
argument  form  of  the  Hankel  functions  given  by 


lim  Ifl2\x)  ->  J^-j^e  ]X 

x—>oc  V  nx 

(4.72) 

which  yields 

lim  hi  \x)  -»  >*/+1  ; 

x— »oo  X 

(4.73) 

therefore, 

(kRhW(kR))1  .ue-JkR 
JiJoo  kR  kR 

(4.74) 

At  this  point,  specializing  the  solutions  to  the  case  of  plane 

wave  incidence  by 

moving  the  point  source  far  away,  it  is  noted  that 

e  ~jkd 

E'  =  >**'  M  ■ 

(4.75) 
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Applying  this  to  Equations  (4.68)  and  (4.69)  yields 

21* 


(4.76) 

. 

=  Sint  (  N 

(4.77) 

where 

i 

j 

i 

•r’K'H 

(4.78) 

! 

u+lV'(v  +  l){2v  +  l)jl,{kRo) 

*  v  o)  *Ro  ' 

(4.79) 

The  expressions  for  the  Ar  and  agree  with  those  given  in  Bowman  et  al.  [33). 

In  Equations  (4.64)-(4.79),  it  is  clear  that  the  sum  is  over  the  eigenvalues 
given  by  solving  Equations  (4.6)  and  (4.8).  All  of  the  previous  expressions  are 
valid  only  for  axial  incidence. 

In  the  next  section,  the  behavior  near  the  tip  and  the  dominant  behavior  in 
the  tip  region  are  determined. 


4.5  FIELD  BEHAVIOR  NEAR  THE  CONE  TIP 


Looking  at  Equations  (4.26),  (4.32),  (4.37),  and  (4.38)  from  the  previous  sec¬ 
tion,  the  fields  could  be  written  simply  as 


E(R)  =  ^MjV)  +  ^N^fR) 

P  <7 


(4.80) 


H(R)  =>K0^apNj,1)(R)  +  ./Yo^T&XV) 


(4.81) 


where  Yq  is  the  characteristic  admittance.  The  superscript  (l)  implies  that  the 
radial  behavior  is  characterized  by  j„(p). 

The  behavior  near  the  tip  is  desired  As  p  ~  kR  — »  0,  the  small  argument 
form  of  the  Bessel  function  is  used.  It  is  giver;  by 


JM 


\ 

2l'+1r (i/  -  3  2) 


OivP 


(4.82) 


Using  this  representation  for  the  Bessel  functions,  and  writing  M  and  N  in  terms 
of  their  auxiliary  wave  functions  yields 


Mp%)  o^/mp  (4.83) 

N^(p)  «  +  1)  x  mq)  (4.84) 

Inserting  these  asymptotic  relationships  into  Equations  (4.80)  and  (4.81)  yields 


E(R)  ss  ^  app^nip  +  ^  bqp'1  *7(7  +  1)  yYqR  +  -R  x  mq  j  (4.85) 

n  n  V  T  / 


and 


H(R)  «  jY0J2d'pPU~1‘/(,/  + 

p 


YpR  +  -R  x  mq 
v 


jy< oX^v7m<i  (4-86) 


where 

cip  —  dpdi/  (4.87) 

bq  =  bqa (4.88) 

We  can  determine  a  bound  on  the  values  of  u  and  7  by  imposing  the  tip 
condition  as  R  — »  e  given  by 


lim  t2  /  Ex  H*  ■  Rd  0  — *■  bounded  at  the  origin 
e— *o  yn 


(4.89) 


where  fi  is  a  shrinking  spherical  volume  centered  at  the  cone  tip.  After  performing 
the  integration,  and  evaluating  the  limit,  the  powers  of  1/  and  7  for  which  the 
integral  would  be  bounded  are  determined.  In  order  to  satisfy  this  condition 


4 


v  >  2 


^ >  2 


(4.90) 


(4.91) 


Closer  investigation  of  Equation  (4.85)  shows  the  dominant  behavior  for  the 


E  field  given  by 


Er-p' 

Ee^p'1 


E&  P1 


(4.92) 

(4.93) 

(4.94) 


where  7  is  given  by  Equation  (4.8). 


The  surface  current  behavior  is  determined  by  working  with  Equation  (4.86). 


The  dominant  behavior  is  found  to  be 


sR  P 


1*4*  -  P* 


(4.95) 


(4.96) 


where  u  is  given  by  Equation  (4.6).  Now,  the  dominant  tip  current  will  be  inves¬ 
tigated  in  more  detail. 

Returning  to  axial  incidence  case,  evaluating  n  x  H  on  the  surface  at  Rq  = 
(/?0,  Qq,  <t>)  and  substituting  Equation  (4.82)  into  Equations  (4.78)  and  (4.79),  yields 


(4.97) 


■  4.*  «_*  -  «  4  .  4.  #  ,  -  « 


r  ■**  *  4  - 


where 


k Rsa 


Y 


j'(>/2)>  +  l)v^(**o) 


1/-1 


D>/2>’ 


sintf0r  (v  +  l/2)Ll{0o) 
s/*(kRor 


=  £ 


rfr  +  1/2)4  (<?0) 


(4.99) 

(4.100) 


I>  +  1/2)4(0O) 

where  the  ”sa”  means  that  the  small  argument  approximation  for  the  Bessel  func¬ 
tion  has  been  applied. 

When  90°  <  0q  <  180°,  0  <  v\  <  1.  Specifically,  when  0q  is  165°,  U\  =  .96714 
and  ix  =  1.03163.  Since  IcRq  <<  1,  it  is  clear  that  i/x  is  the  dominant  eigenvalue. 
Therefore,  the  dominant  term  for  the  current  with  axial  incidence  is  given  by 


i*'  in  \  ~2jT°  <  «»*  '  (j/2)1-1 11'  1  +  1  in.| 


6o 


floHi'i  +  1/2)4,  (^o) 


v  ~2jY°  I  -sin^  1  O/2)1^1  *^1(^1  +  l)v/*(fcfloP  1  (i  lfvM 

J^(Ro)  "  X  cos 4>  i  TUYYU2)Ln0n)  (4-102) 


sin l  cos0  )  r(t/!  +  1/2)4,  (0O) 

This  is  also  true  for  non-axial  incidence,  since  m=l  still  produces  the  eigen¬ 
value  for  which  <  1.  It  must  be  modified  to  account  for  the  fact  that  /?  ^  0. 
The  expressions  are  given  by 


>  =  Y  CM  f  COS  <t>  1  1  xYUkRo)^  1 

Jji;(Ro)  '  ^sin2^  I  sin<£  J  2-i  Hr/,  +  l/2)L,U0n)  (  '  ’ 


U\  snr  0q 


2-1-2  r(^  + 1/2)11,(00) 


J>o)  ■  { 


G{0)  j  -  sin  0 


}-L 

)  2-1- 


y/n{jkR0)^ 


-1 


cos<^  J  2-1-2  r(i/,  +  1/2)4, (0o) 


(4.104) 
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where 


j  ph  (cos  P)  j 

G[0)  =  <  sin/?  >  .  (4.105) 

l  &W)  J 

The  upper  and  lower  terms  in  (4.105)  are  associated  with  the  O'  and  4>'  notation 
from  above.  These  expressions  also  agree  with  the  expressions  for  the  dominant 
behavior  near  the  tip  found  in  Bowman  et  al.  [33].  In  the  next  section,  the  dominant 
tip  current  is  compared  to  the  GO  current  for  the  case  of  axial  incidence. 

4.6  COMPARISON  OF  THE  DOMINANT  TIP  CURRENT  TO  THE 
GO  CURRENT 


The  dominant  tip  current  and  the  GO  current  will  be  compared  for  the  case 
of  axial  incidence.  It  will  be  shown  that  these  currents  are  very  similar  in  this 


case. 


As  mentioned  previously,  from  Equations  (4.103)  and  (4.104)  (replacing  Rq 
with  r'  for  this  discussion)  since  kr'  <<  1,  and  j/j  <  1,  it  is  clear  that  ui  is  the 
dominant  eigenvalue.  The  eigenvalue  is  such  that  V\  <  1  for  cone  angles  in  the 
range  given  by  90°  <  Oq  <  180°;  therefore,  the  dominant  term  for  the  current  is 


given  by 


and 


where 


and 


=  y  (  COS  <f)  1  1  y/njjkr')^  1 

sR{  ’  Vjsin2^  l  sin<*>  /2«'i-2r(i/1  +  l/2)Lj1(«0) 

jl'(  i)  =  Y  G(ft)  f  -  sin <j>  )  1  y 

1  °sin0o  1  cos0  /  2t'i-2l>1  +  1/2)1^  (0O) 

(  (cos/?) 

G(0)  =  {  sin/?  • 

I  , 

Ll(0O)  =  d%jPv(cos9  o)- 


(4.106) 


(4.107) 


(4.108) 


(4.109) 
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Combining  these  terms  yields 


k 


k.v. 


J?V)=  Sjf  } +*'.**{  }]  (4.110) 


where 


Jo  =  Y0 


G{0)  1 


(4.111) 


t/\  sin  Oq  2"i  2  r(i/j  +  1/2)  sin  OqL^  (#o) 

The  normalization  integral,  given  in  Appendix  D,  for  this  set  of  eigenvalues  is  given 
by 


/>o)  = 


-sin  60  . 


2t/  +  1 


Pt(c<*9o)Ll{0o). 


(4.112) 


Therefore,  one  can  write 


ninfl  T  i(ff  )  +  1  ^  *<>) 

smSo-M^o)  - - ivf/ - — 

(cos  0O 


(4.113) 


Substituting  this  into  the  expression  for  the  current  (4.110)  yields 


4\r')  =  JiUkrT^W) 


(4.114) 


where 


fe(d>')  =  sin  6q  < 


f  x(cos2  <f>'  +  i/j  sin2  <f>')  +  y(l  —  vx)  sin  <j>' cos<f>' 

+  z  cot  e0  cos  4>'  O'  incidence 


x(l  -  u{)  sin  <t>'  cos  <f>'  +  y(sin2  <j>'  +  cos2  4>') 

+ z  cot  Oq  sin  <j>'  <i>'  incidence 

(4.115) 

is  the  vector  for  the  exact  current,  and 

1  y/n  plx  (cos  ^o)  1 


J  i  =  -Y0- 


G{0 ) 


V1(2i/1  +  1)  2^1  2  r(i/j  +  1/2)  sin^o  ll{0 o) 

The  GO  current  is  given  by 


(4.116) 


Jgo(r',*')  =  2y0e’*r'‘”%„(<)>') 


(4.117) 
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where 


m 


& 

S 


i 


dS» 


1 

.1 


{x  sin  Oq  z  cos  9q  cos  <//  0'  incidence 

x cos 0q s>n /3 sin 4>'  +  z cos 9q cos 0 sin <f>'  incidence  (4.118) 
+y  (sin  9q  cos  0  -  cos  9q  sin  0  cos  4>') 


cos  7  =  sin  a  sin  0  cos  <f>  +  cos  6q  cos  0 

(4.119) 

For  cone  angles  in  the  range  given  by  90°  <  9q  <  180°,  V\  ss  1. 

If  one  evaluates 

these  expressions  for  axial  incidence  (0  =  0)  and  i/\  «  1,  then 

^,)  f  fshlb  +  iccsW  ^incidence 

(  y  sin  9q  +  z  cos  9q  sin  4>'  4>'  incidence 

(4.120) 

and 

t  ^  |  x  sin  9q  +  z  cos  9q  cos  4>'  91  incidence 

(  y  sin  9q  +  z  cos  9q  sin  <f>'  4>'  incidence. 

(4.121) 

Hence,  one  can  define 

W)  =*«(*')=  fgo(*')- 

(4.122) 

Moreover, 

G(0)=*'l(^l  +  1)  „  J 

(4.123) 

Ji 

r(^i  +  i/2)  ~ 

(4.124) 

PL  (cos  90) 

—  •  ~n -  «  -1 

sin  9q 

(4.125) 

fl(0o)  »  1 

(4.126) 

COS  7  =  COS  9q 

(4.127) 

Lastly,  investigating  the  first  few  terms  in  the  eigenfunction  expansion  with  the 
knowledge  that  for  the  given  range  of  angles 


i/j  «  1 
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(4.128) 


i/2  »  2  (4.129) 

<71  ss  t/2  ~  1  83  1  (4.130) 

llx{* o)  «  1  (4.131) 

/‘2(*o)  *  2  (4.132) 

/-},(*<))  «  1,  (4.133) 

one  can  see  that  the  series  is  that  of  an  exponential  given  by 

ss  C(1  -f  jkr' cos  0qcos /?  +  ...].  (4.134) 

Therefore,  for  axial  incidence,  this  exponential  is 

(4.135) 

This  is  exactly  the  same  phase  as  that  given  for  the  GO  current.  The  two  currents 
may  be  written  as, 

Jgo(r',*')  =  2roejfcr'co8®°f(<*')  (4.136) 

and 

=  ly0(jVp-V*r'c“,0f(*').  (4.137) 

The  [jkr,)v i-1  has  been  kept  in  tact  to  preserve  the  singularity  at  the  tip. 

The  behavior  of  these  two  currents  is  very  similar  for  the  axial  incidence  case; 
therefore,  it  is  felt  that  the  matching  of  the  current  Je  near  the  tip,  to  the  current 
Jgo  far  from  the  tip,  could  provide  a  very  good  approximation  to  the  actual  current. 
It  is  proposed  as  a  part  of  future  research,  that  this  matching  be  found  in  terms 
of  the  transition  function  F  from  (3.30)  with  a  proper  argument.  Such  a  current 
representation  is  expected  to  be  more  accurate  than  the  one  based  on  PO.  Once 
this  more  refined  approximation  for  the  current  is  found,  one  could  calculate  the 
tip  diffracted  field  using  the  asymptotic  methods  in  Chapter  III. 
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In  the  next  section,  reciprocity  is  invoked  to  determine  the  far  fields  of  a 
circumferential  slot  on  the  cone  using  the  Green’s  function  derived  in  this  work. 

4.7  RADIATION  FROM  CIRCUMFERENTIAL  SLOT  ANTENNAS 


The  Lorentz  Reciprocity  Theorem  as  shown  in  [34]  relates  the  fields  from  one 
set  of  sources  to  those  of  another  set  of  sources.  That  theorem  is  used  here  to 
find  the  far  zone  fields  of  a  slot  in  a  cone  in  terms  of  the  reciprocal  problem  of 
determining  the  current  induced  on  the  cone  surface  by  a  distant  electric  current 
point  source;  the  electric  current  point  source  is  located  where  the  far  field  of  the 
slot  is  to  be  evaluated,  and  the  slot  is  replaced  by  an  equivalent  magnetic  surface 
current  distribution  over  the  slot  aperture  with  the  slot  short-circuited.  It  is  noted 
that  the  solution  to  the  reciprocal  problem  of  determining  the  current  induced  on 
the  cone  by  the  distant  electric  point  source  can  be  found  from  the  previous  two 
sections.  In  this  case,  these  fields  radiate  in  the  presence  of  a  perfectly  conducting 
semi-infinite  cone.  Reciprocity  states  that 


1 1  jv  [EaJb-HaMb]dv  =  J  J  jv  [Eb  ■  Ja  —  Hb  •  Ma]  dv  (4.138) 

where  (Ea,Ha)  and  (Eb,Hb)  are  the  fields  radiated  in  the  presence  of  the  cone 
by  sources  (Ja,Ma)  and  (Jb,Mb)  respectively.  This  is  specialized  such  that  we 
have  either  an  electric  or  a  magnetic  source.  Let  (Em,Hm)  be  the  fields  radiated 
by  M,  and  let  (Ee,He)  be  the  fields  radiated  by  J.  Therefore, 


Ja  =  0 

(4.139) 

Jb  =  J 

(4.140) 

Ma  =  M 

(4.141) 

Mb  =  0 

(4.142) 
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In  this  case,  the  reciprocity  integral  yields 


I  j  Em(Rf)  •  J(R/)dVr/  =  ~  J  I  Jv  HelRO-MtR'JdV''  (4.145) 


Defining  the  electric  current  as 


J(R')  =  pe(R)6(R  -  R') 


(4.146) 


then 


Em(R)  •  Pe(R)  =  -  j  I  fv  He(R/)  •  M(R,)dVr/.  (4.147) 

The  Green’s  function  given  previously  is  for  electric  sources.  From  Equation  (4.16), 

Ee(R')  =  jfcZcf(R,,R)  -pe(R). 


Applying  Maxwell’s  equation  yields 

He(R')  =  -V'  x  f  (R',R)  •  pe(R). 
Substituting  into  Equation  (4.147)  yields 

Em(R)  ■  Pe(R)  =  J  j  jv  V'  x  f(R',R)  -pe(R)  •M(R,)dV’/ 
In  [32],  it  is  given  that 

V'  X  f(R',R)  •  Pe(R)  =  Pe(R)  •  V'  X  f(R',R) 

but 

V'  x  f(R',R)  =  V  x  f(R,R'). 

We  can  write  Equation  (4.150)  as 

Pe(R)  ■  Em(R)  =  pe(R)  •  /  /  jv  v  X  f(R,R')  •  M(R ')dV'. 


(4.148) 


(4.149) 


(4.150) 


(4.151) 


(4.152) 


(4.153) 


The  dipole  moment  pe(R)  was  arbitrary;  therefore,  (dropping  the  m  subscript  for 


convenience) 


/ 

E(R)  =  f  f  f  V  xf(R,R')-M(R ')dV'.  (4.154) 

Now  specialize  this  volume  current  M  to  a  surface  current  Mg.  Letting  an  E  field 
exist  only  within  the  slot,  the  voltage  can  be  written  by  V  =  -E-dR.  This  voltage 
only  exists  at  R'  —  a;  therefore, 


E8iot(R')  =  ~RV6(R'  -  a). 


(4.155) 


But  the  equivalent  magnetic  surface  current,  Mg,  is  needed.  This  is  given  by 


Mg(R')  =  -Eglot(R')  x  0  =  V S(R'  -  a)*' 


(4.156) 


and  dS  =  R'  sin  OodR'd<}>'.  Therefore, 


E(R)  =  J  x  T(R,R')  •  V6(R'  -  a)t'R' si 


sin  0q  dR'd<f>. 


(4.157) 


Our  excitation  is  symmetric;  hence,  there  is  no  4>  variation.  This  means  that 
m  =  0.  Therefore, 

E(R)  =  2nVasm0oV  x  f(R,R')  •  ^']^=o ^,=  Q  (4.158) 


The  next  step  is  to  look  at  the  Green’s  function  which  is  given  by 

f  ^  M<V)M(4){R')  ,  _  N<V)N<4)(R')  „  . 


f  ^  M<V)M<4){R')  ,  ^  nS1}(R)n(4)(R')  „  . 

fm  n')=  J  E p  — ^  +  "V  “  R<R 

’  r  M^lRJM^lR')  r  N<4)(R)N<l)(R')  ^  , 

{  LP  - Tip - +  - TT -  R>R 


(4.159) 


We  need  V  x  r(R,R/),  where  V  operates  on  the  unprimed  variables.  Recalling 
that  V  x  M  =  IfcN  and  V  x  M  =  fcN,  yields 

f  ,<W4V>+tr  ^VKV)  R<RI 

T7  V  TVTJ  w't  -  )  nP  ...  9  .  . 


Vxr(R,R')=  F  ,4)  *‘p(1)  ,  *  (4)  9  fl)  , 

1 1 Ep  +  k E? R>i?- 

(4.160) 


L 
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w 
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Next,  performing  the  dot  product  with  <t>'  yields 

Mf"(R')  ■  „,m-  0 =  -Mk^PAtost o) 


N’"‘RMm=()  =  0 


(4.161) 


(4.162) 


Therefore,  there  is  only  a  value  for  p.  It  is  given  by 


p  =  eOu. 


(4.163) 


An  important  thing  to  note  is  that  p  is  determined  by  satisfying  the  boundary 


condition  given  by 


which  yields 


0  x  V  x  r(R,R')  =0 


i  *  N^(R)  =  o. 


(4.164) 


(4.165) 


Therefore,  i/  is  determined  from  the  solutions  to 


P^cosflo)  =  0. 


(4.166) 


Hence,  the  electric  field  is  given  by, 


E(R)  =  2nVkas\n90Y"  — -(-^(fcaJl^rP^costfo)-  (4.167) 

V  ‘‘01/ 

We  are  interested  in  the  far  field.  Therefore,  only  the  transverse  components  are 
needed.  These  are  given  by 

n!;>(R)  •  (7  -  RR)  =  (4.168) 


For  these  boundary  conditions,  the  normalization  factor  is  given  by 
_  .  -  27T  v[v  +  1)  d  _  ,  .  .  d 

=  smhTk^TTdi  p"(  cosS°)S^(cos*°)- 


(4.169) 


1  v  V  V'.-  V  V  V 

vVaav.a, 


Evaluating  the  large  argument  form  of  the  Hankel  function  according  to  the  results 
given  by  Equations  (4.72)-(4.73).  and  recalling  that 


1/2  f1)' 


(4.170) 


and  after  making  the  appropriate  substitutions,  the  expression  becomes 


E„(R)  =  E0^f 


2  v  +  1 
v{v+  1) 


V+l 


/2{ka) 


A^(cosg) 

^P^(cos^o) 


(4.171) 


where 

nr  e~]kR 

E°=vHwa-w  <4172> 

This  expression  matches  that  given  in  both  in  [3]  and  [34].  A  pattern  has  been 
calculated  for  ka  =  n,  Oq  =  165°  which  duplicates  that  given  in  [3],  The  comparison 
of  the  previous  result  with  that  found  using  the  present  work  is  shown  in  Figure 
25.  In  addition,  a  radiation  pattern  has  been  plotted  for  the  case  when  ka  =  10 
and  $ o  =  170°.  This  plot  agrees  with  that  given  in  [35]  and  therefore  supports 
their  contention  that  some  of  the  figures  shown  in  [3]  are  in  error.  The  comparison 
of  the  previous  result  and  that  found  using  the  present  work  is  shown  in  Figure 


Figure  25:  Comparison  of  Radiation  Pattern  for  Circumferential  Slot  at  ka 

and  8q  =  165°  from  [3] 


Figure  26:  Comparison  of  Radiated  Pattern  for  Circumferential  Slot  at  ka  =  10 

and  9q  -  170°  from  |35] 
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CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


This  dissertation  has  endeavored  to  determine  an  expression  for  the  plane 
wave  fields  scattered  by  a  fully  illuminated,  semi-infinite,  perfectly  conducting 
cone  such  that  it  remains  valid  throughout  the  region  exterior  to  the  cone.  To 
this  end,  a  uniform  asymptotic  representation  based  on  the  physical  optics  (PO) 
approximation  was  developed.  Using  these  expressions,  a  UTD  tip  diffracted  field 
was  determined.  With  the  use  of  PO,  however,  comes  the  inherent  shortfalls. 
Namely,  the  fact  that  PO  does  not  satisfy  the  proper  boundary  conditions,  and  that 
it  does  not  properly  model  the  tip  discontinuity.  In  spite  of  these  shortcomings, 
the  results  are  very  accurate  as  seen  in  comparison  with  the  rigorous  asymptotic 
solution  for  the  narrow  angle  cone  given  in  [5].  Although  the  present  analysis 
was  performed  for  the  semi-infinite  cone,  it  is  also  valid  for  the  finite  cone  since 
high  frequency  diffraction  is  a  local  phenomenon.  Recall  that  the  previous  GTD 
solution  for  the  RCS  of  the  flat-backed  cone  [10]  had  not  accounted  for  the  tip- 
base  interaction.  To  verify  the  validity  of  the  UTD  tip  diffracted  field  given  in  the 
present  work,  it  is  shown  that  when  it  is  included  in  the  previous  GTD  solution  for 
the  RCS  of  a  finite  flat-backed  cone  [10],  then  it  serves  to  correct  the  long-standing 
discrepancy  between  calculation  and  measurement.  Although  the  application  of 
the  UTD  tip  diffracted  field  was  concentrated  on  the  scattering  problem,  it  is 
intimately  related  to  the  radiation  problem  through  reciprocity.  It  could  be  used 


v  -J- 


fT»  U  '  V  f  7« 


equally  well  to  determine  the  fields  radiated  by  an  antenna  on  the  cone  surface. 

For  the  case  of  axial  incidence,  the  actual  tip  current  and  the  GO  based  current 
were  compared  and  it  was  shown  that  the  two  currents  are  very  similar  in  this  case. 
Since  the  dominant  term  from  the  eigenfunction  expansion  for  the  current  exhibits 
the  proper  singularity  at  the  tip,  and  the  GO  current  dominates  far  from  the  tip, 
it  is  felt  that  the  matching  of  these  two  currents  in  their  respective  regions  by  a 
transition  function  with  a  proper  agrument  could  provide  a  better  model  for  the 
actual  current.  The  fields  from  this  current  could  be  evaluated  using  the  method 
given  in  Chapter  III. 

The  necessary  computer  routines  for  determining  the  eigenvalues  associated 
with  the  exact  eigenfunction  solution  for  the  cone  problem,  as  well  as  the  computer 
routines  needed  to  calculate  the  associated  Legendre  polynomials  evaluated  at 
these  eigenvalues  have  been  written  and  tested  by  computing  the  far  fields  radiated 
by  a  circumferential  slot  on  the  cone.  In  performing  this  check,  the  results  were 
used  to  verify  that  some  of  the  plots  previously  publisned  in  [3]  are  in  error. 

There  are  several  recommendations  for  future  research.  Now  that  the  form  of 
the  tip  surface  wave  is  better  understood,  a  hybrid  MM/GTD  analysis  could  be 
undertaken  to  further  refine  the  tip  diffraction  mechanism  along  the  cone  surface. 
Also,  ways  to  modify  the  PO  solution  to  enforce  the  boundary  conditions  could 
be  investigated.  The  case  when  the  cone  is  not  fully  illuminated  needs  to  be 
investigated  also.  In  this  case,  the  creeping  wave  in  the  shadowed  region  of  the 
surface  needs  to  be  accounted  for,  and  the  range  of  PO  integration  which  is  dictated 
by  the  limits  of  the  lit  region  on  the  surface  must  be  modified  via  Fock  theory.  In 
addition,  the  PO  correction  based  on  the  dominant  term  of  the  exact  eigenfunction 
could  be  tested  and  improved  upon.  The  availability  of  the  computer  routines 
required  to  calculate  the  exact  eigenfunction  solution  enlarges  the  possibilities  for 
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research  that  integrates  this  eigenfunction  solution  with  high  frequency  techniques. 
Last,  but  not  least,  the  asymptotic  evaluation  over  the  azimuth  variable  (</>')  in  the 
development  of  a  UTD  solution  as  in  Chapter  III  should  be  investigated.  This  was 
investigated  in  a  cursory  manner;  however,  there  were  some  problems,  and  due  to 
time  constraints,  not  completed.  One  should  be  able  to  determine  the  stationary 
phase  contribution  that  is  the  source  of  the  reflected  field. 


APPENDIX  A 


EVALUATION  OF  TERMS  IN  THE  UNIFORM  SADDLE  POINT 

INTEGRATION 

A.l  INTRODUCTION 

In  this  appendix,  the  terms  required  to  evaluate  the  uniform  saddle  point 
integral  are  shown. 

A. 2  DERIVATION  AND  EVALUATION  OF  TERMS 

Transforming  to  the  steepest  descent  path  required  that 

CM  =  f(r',*')^.  (A.l) 

This  must  be  evaluated  at  s=0  and  s  =  sa.  For  s=0 

G(0)=f(r.„«')^|r,=r<,  (,1.2) 

and  for  s  =  sa 

CM)  =f(0,/)^  (A3) 

but 

f(0,<fi')  =  0.  {A.  4) 

Hence, 


G(Sa)  =  0. 


(A.S) 


I 

! 


I 


Turning  to  H(s),  recall  that 


H(s)  =  1  [<?(«)- 0(0)1  =  »c'(«)  -  0(3)  +  0(0) 


Evaluating  this  at  s  =  0,  yields  Applying  L’HospitaPs  rule  yields 

zr.ni  _  G"(0) 


H(  0)  = 


When  evaluated  at  s  =  sa, 


H{sa)  =  ^g,(^a)+g(0) 


(A.6) 


(A.7) 


(A.S) 


In  general, 


j2  _/  j  j 

G'W  -  t(r',4,')—+ty^)^  (A.9) 

G"w  =  (*•“) 

The  next  task  is  to  determine  the  expressions  for  the  differentials  above.  Recall 


that  the  transformation  was  given  by 


q{r',tf>')  =  q{rs,<j>)  -  s2 


(A.ll) 


This  is  used  to  derive  the  expressions  for  the  differentials  in  terms  of  q  and  its 


Recall  that  q'^g,#')  =  0,  evaluating  these  for  r'  =  r3,  yields 


—  ] 
ds  lr=r* 

dV, 


ds2  lr=r* 


V  «*('.) 

2  q'"(ra) 
3[9"(ra)l2 


dv  1  -  1  1 

/  -2 

Li.  5[«'"(r,)|2l 

ds3  lr=r‘  2 [g"(ra)  2  V 

9"(r*) 

[q  3  q»(rs) 

(A. 16) 
(A.17) 
(A.18) 


I  -2 

l"(rs) 


(A. 19) 


These  are  then  used  in  the  expressions  for  G(0)  and  H(0).  These  are  now  given  by 

GT(0)  =  f(r.)1 

and 

Turning  to  those  evaluated  at  sa ,  one  obtains 


5 

f?W(rs)l 

2  r(r,)] 

3 

l  *"('.) J 

«"(r.)  J 

(^.20) 

dr 1  .  -2  sa 

dl'r'=°~  7(0)' 


(A.2I) 


It  has  already  been  shown  that  G(s0)  =  0,  Next,  G'(sa)  is  determined.  Recall 


and  since  f(0)  =  0, 


Lastly, 


(A.  22) 

0, 

W/  \  4^(°)5o2 

(A. 23) 

G"M  =  f"{0,7)(d/  +  3t'(0,7)ds2ds 

(A. 24) 

where  the  differentials  are  also  evaluated  at  r'  =  0.  Therefore,  expressions  are 
needed  for  this  case.  One  should  note  that  if  the  function  q  is  expanded  in  a 
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0,  the  differentials  for  the 


Taylor  series  about  rs,  it  can  be  shown  that  as  rs  — + 
r'  =  0  case,  which  are  shown  below,  reduce  to  those  given  for  r1  =  ra. 


—  1  /  - 

2  sa 

ds  =°  - 

^(0) 

d2r>  \ , 

2 

ds2  Jr -°  ~ 

k'(o)] 

rfv  1 ,  - 

4 

ds3  M=0  ~ 

[,'(0)]3 

{2«aV'(0)  +  |?'(O)|2} 

r?"(o)f 


2s 


—  6sa 


-  3sa9w(0) 


(A.25) 
(A.26) 
(A. 27) 


a  ?'(o)  —  L?'(o)  J 

Expressions  for  the  phase  functions  used  in  the  differentials  are  needed.  More¬ 
over,  the  expressions  for  f,  f',  and  f"  are  needed.  Beginning  with  the  phase  func¬ 
tions.  These  are  given  by 


q(r')  =  r'  cos  7  -  \J r2  r2  -  2rr'  cos  f 

r'  —  r  cos  f 


q'(r' )  =  COS7- 

«V)  = 


\J r2  +  r2  —  2  rr1  cos  f 
— r2(l  —  cos2  f) 


«*V)  = 


?‘V)  = 


(r2  +  r 2  -  2rr'cosf)3/2 
3r2(r'  -  rcosf)(l  -  cos2  f) 


(r2  +  r2  -  2rr,cosf)5/2 
3r2  ( 1  —  cos2  f) 

(r2  +  r2  -  2  rr1  cos  f)7/2 


r2(l  —  5 cos2  f)  +  8r/cos  f  —  4 r  2 


(A. 28) 
(A. 29) 

(A. 30) 
(A.31) 

(A. 32) 


Evaluating  these  at  r1  =  0  and  r 1  =  rs  yields 


9(0)  =  -r 

(A. 33) 

q[rs)  =  rcos(f  +  7) 

(A.34) 

q'( 0)  =  cos  f  + cos  7 

(A. 35) 

9'(ra)  =  0 

(A. 36) 

,"(0)  = 

(A.37) 
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r  sine 


(A.38) 


<"'(<>)  =  -3 


cos  e  sin  e 


(A.39) 


m,  ,  3  cos  7  sin4  7 

9  (r‘)  =  12 - TTT- 

rL  sin  e 


(A.40) 


-3  sin5  7  2 


9*U(r*)  =  TT-3-{4cos  7~sinz7] 

rd  sm^e 


(A.41) 


These  could  also  be  found  by  noting  that  when  ra  =  0,  then  f +  7  =  n.  Substituting 


this  into  the  equations  for  r'  =  r3  the  values  reduce  to  those  shown  for  r  =  0. 


To  find  G(0)  and  H(0)  we  need 


j"(rs)  0  sin  7 
q»(rs)  7  r  sine 


(.4.42) 


9M(ra)  r2sin‘,e 

Therefore  G(0)  and  H(0)  are  given  by 


(A. 43) 


,  /2r  sin  e 
G{0)=t{ra)J—r± 
V  sin  7 


(A.44) 


H(0)  =  i  (^)S/2  Ur.)  +  +  J  .  (A45, 

2  \  sin  7  /  rsinf  4rzsm  e 


The  remaining  task  is  the  derivation  of  the  expressions  for  f  and  its  derivatives. 


Recall  that 


f  {r1  ,<f>')  =  +  — )k  x  (xsin  a  -  zcosacos<j>'). 

2  7r  R  R 


(A.46) 


Since 


R-* 
R  ~  R' 


(A.47) 


wmmtmm  mmm: 


we  write 


f(r\*')  =  g(r')  x  h(A 


(A- 48) 


where 


g(r')  =  P{r')R, 

(4.49) 

with 

,  ,ybsina,  ,  1 , 

,,(r)  =  -r  {3k  +  Rh 

(>1.50) 

and 

h  {4>)  =  xsina  —  z  cos  a  cos  <f> . 

(>1.51) 

Taking  derivatives  of  f  requires 

f(r'J)  =  ^  (g(r')  x  h(*')] 

(>1.52) 

and 

f'W)  =  [s(r')  X  h(*')l  . 

(4.53) 

Since  h(0/) 

is  not  a  function  of  r1, 

f W)  =  *V) *  MA 

(4.54) 

Similarly, 

f"(r',^')  =  g"(r')  x  h(</>'). 

(4.55) 

Recall  that 

g(r')  =  P(r')R. 

(4.56) 

Therefore, 

*V)  =  p'(r')R  +  p(r')^ 

(4.57) 

and 

g'V)=P'V)R  +  2pV)^+p(r')~, 

(4.58) 

(>4.58) 


where 


(4.59) 


(>1.60) 

(A.61) 

(4.62) 

(>1.63) 

(A.64) 

(4.65) 

(4.66) 

(4.67) 

(4.68) 

(4.69) 

(4.70) 

(4.71) 


and 


V2  =  f  ^sin2  a  cos  0  sin  <f>'  +  sin  a  cos  a  sin  6  (sin  <f>  —  cos  4>  sin(^  — 

$  |  —  sin2  a  sin  0  sin  <t>  +  sin  a  cos  a  cos  6  (sin  <f>  -  cos  ^'sin(^  —  ^))J 
4>  sin  a  cos  a  [cos  4>  -  cos  4>  cos [4>  -  <f>')]  (A. 72) 

Now  the  expressions  for  p  and  its  derivatives  are  shown.  These  will  be  evaluated 


at  the  saddle  point  r'  =  r4,  and  the  end  point  r'  =  0. 


P(  0)  =  0 

,  ,  „,sina  sin2  7  sin  7 

pr*  =  -Y0——ra  —  f-  jk+— r-L 
2n  r2  sin^f  rsmfj 

sina  ll 

p(0)  = 

,'M  -  final 

2n  r2  sinz?  IL  rsinfj 

sin 7  r  ..  „  sin 7  ) 

-rscos7 — ; —  \2jk  +  3 — : —  > 
rsinf  L  rsm  ?  J 

sina  T  3l 

P  (0)  =  -Yq — =-  \2jk  +  -  cos 5 


(A. 73) 
(A.74) 


(A. 75) 


(A. 76) 


(A. 77) 


\  ,,  sina  sin3  7  (  L.,  „  sin^  1  f  sir 

p  (r3)  =  -y0— - r  ;-o~  {  -  2jfc  +  3— 7—  2cos7  +  r,— : 

2n  r3  sin3  f  (.  rsinf  rsi 


sin  7  .0 
— —  sin^  7 
’sin  f 


sin^  2  sin  7  ) 
+3rs  2jk  +  4 — ; —  cos^7 — 7 — > 
r  sin  f  r  sin  f  J 


(A.78) 


As  with  the  phase  functions,  when  rs  =  0,  then  f  +  7  =  *.  Under  this  condition, 
the  expressions  evaluated  at  r'  =  rs  yield  the  expressions  for  those  evaluated  at 


r  =  0. 


The  expressions  in  this  appendix  were  used  to  program  the  results  of  Chapter 
III  for  the  uniform  PO  solution. 


APPENDIX  B 


VECTOR  WAVE  FUNCTIONS 

B.l  INTRODUCTION 

In  this  appendix  the  derivation  for  the  vector  wave  functions  in  the  general 
case  is  shown.  This  is  then  specialized  to  the  conical  geometry. 

B.2  GENERAL  CASE 

In  a  source-free  region  of  space,  both  E  and  H  fields  satisfy  the  vector 
Helmholtz  equation  given  by 

Vx(Vx{h  })-*2{h  }  =  °  (Bl) 

As  shown  in  [8,20,21],  the  vector  wave  functions  are  also  solutions  to  the  vector 
Helmholtz  equation  and  are  defined  by 

M9  =  V  x  ipqf  =  x  N9  (5. 2) 

N,  =  iv  x  M9  (5.3) 

where  is  the  solution  to  the  scalar  wave  equation 

(V2  +  k2)xl>q  =  0  (5.4) 

and  f  is  a  vector  function  to  be  determined.  Substituting  Equation  (B.2)  into 
Equation  (B.l)  yields 


By  definition,  V  x  Vu  =  0;  therefore,  we  define 

V  x  V  x  4>f  —  k2ipi  =  Vu  (5. 6) 

Using  a  vector  identity,  yields 

V  x  (V  x  ipf)  =  V  x  (Vtp  x  f  +  ipV  x  f)  (J5.7) 

Looking  at  each  term  separately,  it  can  be  seen  that 

V  x  (V0  x  f)  =  (f  ■  V)V0  -  V20f  -  V0  •  Vf  +  V^(V  •  f)  (B.8) 

and 

V  x  (ipV  x  f)  =  Vip  x  (V  x  f)  +  0(V  x  V  x  f)  (5.9) 

After  expanding  Equation  (B.9),  adding  it  to  Equation  (B.8),  and  performing  a 
fair  amount  of  vector  algebra 

V  x  (V  x  xl>{)  =  V(V  •  rPf)  -  fV2tl>  -  rl>V2 f  -  2 (V0  •  V)f  (5.10) 

Manipulating  Equation  (B.6)  a  little  further,  yields 

-  V>V2f  -  2(V  •  V)f  -  f(V20  +  k2xjj)  =  V(u  -  V  •  0f)  (5.11) 

Recall  that  ip  satisfies  the  scalar  wave  equation  given  in  Equation  (B.4);  therefore, 

V>V2f  +  2{Vip  •  V)f  =  VU  (5.12) 

where 

U  =  V  ■  ipf -u  (5.13) 

is  another  scalar  function. 

In  order  to  arbitrarily  satisfy  the  boundary  conditions  for  different  ip,  ip  and 
f  must  be  independent.  Since  ip  and  f  must  be  independent,  closer  examination  of 
Equation  (B.12)  yields  two  equations 


2(V0- V)f  =  VC/ 


Equation  (B.15)  can  be  expanded  for  the  x  component  of  f  as 


drpdfx  d\l>dfx  drpdfx  1  dU 


dx  ax  ay  dy  dz  dz  2  dx 


The  expansions  for  fy  and  fz  follow  similarly.  For  independence  we  need  fx  to  be 
only  a  function  of  x,  fy  to  be  only  a  function  of  y,  and  fz  to  be  only  a  function  of 
z.  Also  needed  is  U  =  2 crp.  Enforcing  the  preceeding  conditions  leaves 


The  solutions  to  these  equations  yields 


fx  =  cx  +  al 


fy  =  cy  +  a2 


fz  =  cz  +  a3 


Combined  into  vector  notation,  this  yields 


f  =  c(xx  +  yy  +  zz)  +  a\x  +  a2y  +  a3z 


which  can  be  written 


f  =  cR  +  aa 


It  can  be  shown  that 


V  x  f  =  0 


(£.17) 


(£.18) 


(£.19) 


(£.20) 


(£.21) 


(£.22) 


(£.23) 


(£.24) 


(£.25) 


mmi 


V  •  f  =  3c 


(B.26) 


(V0-V)f  =  cV0 


(B.  27) 


These  relationships,  along  with  Equations  (B.4)  and  (B.24)  yield  the  general  form 
of  the  vector  wave  functions  as 

M  =  -f  x  Vtp  {B. 28) 

N  =  krpf  +  ^V(c^>  +  f  •  V^)  (B.29) 

In  the  next  section,  the  general  result  is  specialized  to  spherical  coordinates 


and  the  spherical  vector  wave  functions  are  used  to  form  the  conical  vector  wave 
functions. 

B.3  CONICAL  VECTOR  WAVE  FUNCTIONS 

The  cone  can  be  described  in  the  spherical  coordinate  system;  therefore,  we 
specialize  the  vector  wave  functions  to  spherical  coordinates.  For  this  case  f  =  R; 
therefore,  c  =  1  and  a  =  0.  The  spherical  vector  wave  functions  form  the  basis  for 
the  conical  vector  wave  functions  [32].  There  are  a  few  differences.  One  is  that  the 
degree  of  the  Legendre  polynomials  is  fractional  and  determined  by  the  boundary 
conditions  of  the  cone.  Another  is  the  orthogonality  relationships  of  the  auxiliary 
vector  wave  functions.  The  vector  wave  functions  for  this  case  have  the  form 

Mq  =  -R  x  Vtfig  (B. 30) 

Nq  =  fcVR  +  +  R  ■  VrPq)  (B.31) 

In  (22) ,  it  is  shown  that  the  tpq  that  satisfies  the  scalar  wave  equation  in 
spherical  coordinates  when  the  origin  is  included  is  given  by 

4°(o)  =  (B.  32) 


_  t 


where  p  =  kR  and  q  =  5  rrw  . 

The  Yq(6,  <f>)  are  called  the  tesseral  or  sectoral  harmonics  and  are  defined  by 


n(M)  =  pri'osn)  { ^ } 


(B.  33) 


where  the  P™(cosO)  are  the  associated  Legendre  polynomials  of  degree  t/,  and 

order  m.  The  zl^(p)  are  given  by 

spherical  Bessel  function 
spherical  Neuman  function 


Zv  ^  ^  h\}\p)  spherical  Hankel  function,  1st  kind 

"  (2) 

hi  ’(p)  spherical  Hankel  function,  2nd  kind 


JAP 

VAP 


1=1 

1=2 

1=3 

1=4 


(B.34) 


The  e  and  o  refer  to  even  or  odd  behavior  in  the  <j>  variation  described  by  the 
bracketed  terms.  The  e  goes  with  the  top  term  and  the  o  with  the  bottom. 

Substituting  Equations  (B.32)  and  (B.33)  into  Equations  (B.30)  and  (B.31) 
and  performing  a  wealth  of  manipulations,  we  finally  have  the  form  of  the  conical 
vector  wave  functions  given  by 


mJV)  =  4V  )m,(M) 


"W 


where  p  =  kR,  q  =  %  1711/1  and 


lq{0,<t>)  =  v(v  +  l)Yq{e,<t>)R 


(B.35) 
(B.  36) 

(B.37) 


-MM) = { z zt } *  ■ -  **•<*«•)  { zzi }  *  (*»> 

The  19,  mq,  and  R  x  m9  are  called  the  auxiliary  vector  wave  functions.  These 
functions  exhibit  surface  orthogonality  over  a  spherical  sector.  That  is,  for  the 
spherical  sector  defined  by  dfl  =  sin  0d0d<t>  and  fl  =  j  %.  }  we  ^ave  *nte” 

grals  of  the  form 


+ 


} ]  sin 9dBd<t> 


(B.39) 


Due  to  the  orthogonality  of  the  Legendre  polynomials  (shown  in  Appendix  D), 
there  is  only  a  value  when  v  =  u' .  This  yields 


7z /„  °  { l£rtpr(c°s')]2  +  [lp"<cos<,>] } sinM#  =  ^ + 


(£.  40) 


where 


and 


sin  6d0 


(£.41) 


£"»  =  {  2  m  ^  0  (#-42) 

The  relation  used  in  (B.40)  was  found  in  [32].  Using  the  integrals  above,  the 
following  orthogonality  relationships  hold  for  the  auxiliary  conical  wave  functions: 


J  Jn  mpi  ■  mpdU  = 
J  J  {R  x  mp/)  ■  ( R  x  mp)dn  = 


{ 


0 

2ir 

(m 


P/  p' 

u{u  +  1  )Jm„  p  =  p' 


and 


■mi/  P  —  P 


(B.43) 


(£.44) 


Moreover,  these  functions  are  vectorially  orthogonal;  therefore,  they  are  mu¬ 
tually  orthogonal.  That  is,  any  integral  of  two  distinct  auxiliary  vector  wave 
functions  over  a  spherical  sector  is  zero. 

As  described  in  [21,23],  the  vector  wave  functions,  as  shown,  are  complete 
with  respect  to  E  and  H  in  a  source-free  region.  If  a  volume  source  density  were 
to  be  included,  the  Lq  functions  would  be  needed  to  be  included  [21];  however, 
these  will  not  be  discussed  in  this  paper. 

In  general, 


E  —  ^  ^  apMp  +  ^  *  6<jN^ 

p  <7 


(£.45) 


Using  V  x  E  =  — jw/iH  along  with  Equations  (B.2)  and  (B.3)  yields 

H  =  jYq  Y,  «pnp  +  JYo  Y,  69M9  (*•*«) 

P  9 

From  these,  the  representations  for  the  surface  currents  earn  be  found  as  well. 

In  the  next  appendix,  the  dyadic  Green’s  function  for  conical  systems  will  be 
derived  in  terms  of  these  vector  wave  functions. 


APPENDIX  C 


DYADIC  GREEN’S  FUNCTION 

C.l  INTRODUCTION 

In  this  appendix  the  dyadic  Green’s  function  is  derived  and  shown  how  it  is 
applied  in  deriving  the  solution. 

C.2  DERIVATION  OF  DYADIC  GREEN’S  FUNCTION 

The  most  general  way  to  relate  vector  fields  to  vector  sources  is  with  the 
dyadic  Green’s  function.  The  dyadic  Green’s  function  satisfies  the  inhomogenous 
vector  differential  equation  shown  below. 

V  x  V  x  f(R,R')  -  Jfe2F(R,R')  =  -7£(R  -  R')  (C.l) 

where  7  is  the  identity  dyad.  Letting  r(R,R')  •  a  =  G(R,R')  to  simplify  the 
derivation  of  Equation  (C.l),  we  have  the  following  equation, 

V  x  V  x  G(R,  R')  -  fc2G(R,  R')  =  -a6(R  -  R')  (C.2) 

Recall,  the  vector  wave  functions  are  orthogonal  functions  that  form  a  com¬ 
plete  set  over  the  solutions  to  source-  free  vector  wave  equation.  As  in  [24],  surfaces 
as  shown  in  Figure  27  are  defined. 

For  R  >  R',  outgoing  waves  that  will  satisfy  the  radiation  condition  on  £,  the 
surface  at  oo  are  needed.  Therefore, 
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Figure  27:  Geometry  used  in  Derivation  of  Vector  Wave  Function 


G“(R,R')  =  5>“(R')m|;"(R)  + 

P  7 

R  >  R'  (C.3) 

where  the  (4)  implies  that  the  radial  behavior  is  described  by  the  spherical  Hankel 
function.  This  function  satisfies  the  radiation  condition  for  the  time  convention 

Green’s  second  identity  for  vectors  is  given  via  the  Divergence  Theorem  as 

/V{A  •  (V  x  V  x  B)  -  B  •  (V  x  (V  x  A)}dV 
=  /5+E{B  x  (V  x  A)  -  A  x  (V  x  B)}  ndS  (C.4) 

As  discussed,  functions  will  be  chosen  such  they  that  satisfy  the  radiation 
condition  on  E;  that  is,  the  integral  over  the  surface  E  at  oo  will  go  to  zero. 
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m 

i 
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Substituting  B  =  NT,  and  A  =  G ^  into  Equation  (C.4)  with  RH  as  a  dummy 
variable  yields 

/v{GI/n(R",R')  •  (V"  x  V"  x  M^fR")) 

— M^(RW)  •  (V"  x  V"  x  GI/n(Rw,R')}dV" 

=  /sfM^R'')  x  V"  x  Gn(R",R') 


-Gn(R",R')  x  V"  x  M^R")}  •  RdS" 


(C.5) 


Note,  in  the  volume  integral,  both  G1  and  G^  appear.  The  volume  integrand  can 
be  simplified  since  M^(R")  satisfies  the  vector  Helmholtz  equation.  Therefore, 
the  left  hand  side  becomes 

-/  M(^(RW)  •  {V"  x  V"  x  GI/II(R",R')  -  Jk2GI/n(RM,R,)}dV"  (C.6) 

Jv  p  ' 

but  by  Equation  (C.3),  the  bracketed  term  can  be  replaced  by  a6(RM  —  R')  and 

/  M^(RW)  •  a6(Rw  —  R;)dVM  =  M^(R/)  •  a  (C.7) 

by  the  sifting  property  of  the  delta  function.  Since  R  <  R’  includes  the  origin,  we 
chose  M^R'  )  which  is  finite  there  since  the  (1)  implies  that  the  radial  behavior  is 
described  by  a  spherical  Bessel  function  which  is  bounded  at  the  origin.  Equation 


(C.5)  can  be  rewritten  as 


JsfM^R")  x  V"  x  Gn(R")  -  Gn(R")  x  V"  x  M^'(R")}  •  RdS" 


(C.8) 


Substituting  the  series  representation  of  G**  from  (C.3),  and  recalling  the  rela¬ 
tionship  between  M  and  N  given  in  Equations  (B.2)  and  (B.3)  where  Vn  operates 
only  on  double  primed  variables,  yields 

^Ep/s^WiM^R")  X  nJ,4)(R")  -n44)(R")  x  N^R")}  •  RdS" 


+*E,/s«?(R'){m'!)(R")  X  m'4>(R")  -  nJ4|(R")  X  n£>(R")>  ■  RdS” 


=  Mjj'(R')  •  a 


(C.9) 


After  expressing  M  and  N  functions  in  terms  of  their  auxiliary  wave  functions, 


and  performing  the  indicated  dot  and  cross  products,  applying  the  orthogonality 


relationships  what  remains  is 


*£  /  (C.10) 

„JS  P  P 


The  terms  within  the  bracket  that  involve  the  spherical  Hankel  and  Bessel  functions 


can  be  rewritten  as 


-  a!?Wm  =  ivu-w, *?’((>))  =  -4  (c-11) 


where  W  is  the  Wronskian  of  those  two  functions.  Recall  that  p  =  kR,  therefore 


Equation  (C.10)  becomes 


i  "V  •  mPdn"  =  mJ|>(R')  •  a  (c.  12) 


For  p  =  o  this  yields 


(C.13) 


From  this,  aJ«(R')  can  be  solved  for  and  given  by 


«8,(R,)-s^»n.»4V)-» 


27Tt/(t'  +  l)/m„ 


(C.14) 


To  find  the  expression  for  aj^R'),  we  would  proceed  in  a  similar  manner  using 


N^(RW)  in  Equation  (C.5)  with  q  =  $  ™7-  This  yields 


_  N  (T*  *  \  _  jktfn  tvt  ( ^ )  /x%t\  » 

a"(R|-2^h  +  l)/m1N2  (R)'“ 


(c.15) 


Substituting  this  into  Equation  (C.3),  we  get  the  representation  for  G**  from  which 
we  can  easily  see  the  form  of  T  as 

Fm  n'l  rWp'WMf'M  ,  ^ n'Vjn'V') 

T(R,R )  =  ^  fT  +  2-  Y, 

p  p  q  q 

R>  R!  (C.16) 

Moreover,  due  to  symmetry  of  the  Green’s  function, 

S,„  ^M^tRJM^IR')  y-N^'lRlN^lR') 

r(R,R )  =  ^  -  +  2^  T 

p  P  q  q 

R<R'  (C.17) 

where 

Op  =  —  l/(l/  +  l)Imu  (C.18) 

J  fcCm 

2ir 

T,  =  777- lb  +  l)/m,  (C.19) 

and  Ima  is  defined  by  Equation  (B.41)  in  Appendix  B.  Note,  the  notation  above 
is  defined  in  the  following  way, 

m^rjmJ/V)  =  m[L(R)m[L(R')  +  m!L(R)m!L(R')  [C.  20) 

and  similarly  for  N^(R)N^(R'). 

Equations  (C.16)  and  (C.17)  combine  to  define  the  dyadic  Green’s  function  for 
all  R.  The  indexes  p  and  q  are  assumed  to  be  distinct.  This  is  required  to  satisfy 
the  boundary  conditions  on  the  cone.  The  Green’s  function  was  chosen  such  that 


»'  x  r(R',R)]<=9o  =  -»  x  r(R',R)]>=#o  =  0  (C.21) 

Since  p  and  q  are  distinct, 

»  x  Mp  (R)]ff_^  =  6  x  N<>(R)]^,o  =  0  (C.22) 
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where  p  =  $  mi/  and  q  =  %  nr y.  Evaluating  these,  it  can  be  seen  that 

6  x  Mp  (R)|  =0  (C. 23) 

J  90 

requires 

«■#)]  =0  (C.24) 

dS  J  9=6 o 

and  that 

hN'1,R)]i.i0=0  (<?•») 

requires 

P7m(cosfl0)  =  0  (C.26) 

The  r  and  7  are  the  theta  eigenvalues.  Their  determination  is  an  important 
part  of  the  problem.  Appendix  D  discusses  this  problem. 

In  the  next  section,  the  procedure  for  using  dyadic  Green’s  functions  to  de¬ 
termine  the  solution  is  discussed. 

C.3  SOLUTION  WITH  DYADIC  GREEN’S  FUNCTIONS 

The  next  step  is  to  show  how  the  dyadic  Green’s  function  is  used  in  the 
solution.  As  in  [25],  beginning  with  the  equation 

V'  x  V'  x  E(R')  -  fc2E(R')  =  -MuIv(R')  (C.27) 

We  introduce  T  that  satisfies  the  primed  version  of  equation  (C.I)  and  is  given  by 

V'  x  V'  x  f(R',R)  -  k2 f  (R',R)  =  -a£(R'  -  R)  (C.28) 

Green’s  second  identity  for  tensors  is  then  applied.  This  yields 

JV{E(R')  •  (V'  x  V'  x  f(R',R))  -  (V'  x  V'  x  E(R'))  -f(R',R)}dV" 

=  /s+eK  x  (V'  x  E(R'))  -f(R',R )dS' 

+  /5+E(n'  x  E(R'))  •  (V'  x  F(R',R))}dS'  (C.29) 


Solving  for  7'  x  V'  x  r(R\R)  using  Equation  (C.28)  and  using  Equation  (C.27), 
to  solve  for  x  V'  x  E(R')  and  substituting  into  the  volume  integral  yields 

-  j  E(R')6(R'  -  R  )dV'  +  juiiJv  Jy{R')  •  F(R',R)dU'  (C.30) 

Also  used  are  the  following  relationships 


V'  x  E(R')  =  — ju>mH(R')  (C.31) 

n'  x  H(R')  =  J8(R')  (C.32) 

n1  x  E(R')  =  -K8(R')  (C.33) 

The  integrals  over  E  vanish  since  E(R')  and  r(R',R)  satisfy  the  radiation  condi¬ 
tion.  Using  the  sifting  property  of  the  delta  function,  and  the  above  relationships 
Equation  (C.29)  can  be  written  as 


+ju>n 


R  i.oU„V  T'  H-^-MR'l  rlR'.RWV" 
J.(R')  •  f(R',R)<iS'  +  jf  K»(R')  ■  (V'  x  F(R',R))rfS' 


(C.34) 


The  form  of  the  equation  can  be  modified  depending  upon  the  material  boundary 
conditions,  and  the  symmetry  properties  of  the  Green’s  function.  This  is  described 
in  many  places,  particularly  in  [21,25,26].  The  Green’s  function  satisfies  the  fol¬ 
lowing  symmetry  relation 

f(R',R)  =  f(R,R')  (C.35) 


In  some  applications  the  volume  integral  is  treated  as,  E*,  the  incident  field.  For 
that  case,  the  expression  for  the  scattered  field,  EB,  is  given  in  terms  of  the  induced 
surface  currents  JS(R)  and  K8(R). 

In  the  next  appendix,  the  derivations  of  the  associated  Legendre  polynomial 


expressions  used  to  program  these  functions  are  discussed. 


APPENDIX  D 


ASSOCIATED  LEGENDRE  POLYNOMIALS 

D.l  INTRODUCTION 

In  order  to  compute  the  exact  solution  to  the  cone  problem  using  the  tesseral 
or  sectoral  harmonics,  we  need  to  be  able  to  calculate  the  associated  Legendre 
polynomials  P™(x)  for  non-integer  degree  u.  Moreover,  in  order  to  determine  the 
required  degrees,  the  solution  to  the  following  are  required 

P™(  cos0o)  =  O  (D.l) 

cos*0)=0  (D. 2) 

with  0o  =  7r  —  a;  where  a  is  the  cone  half-angle.  These  are  solved  as  functions  of  <7 
and  1/  respectively.  For  the  case  of  the  cone,  these  are  generally  non-integer.  This 
is  why  a  special  effort  was  required  to  determine  a  method  to  evaluate  the  Legendre 
polynomials.  In  this  appendix,  T  represents  the  well  known  gamma  function.  We 
will  begin  by  disscussing  the  integral  representation  for  the  Legendre  polynomial. 

D.2  INTEGRAL  REPRESENTATION 

The  solution  of  boundary  value  problems  in  the  spherical  system  gives  rise  to 


the  following  differential  equation 


for  which  one  solution  is  P™( x),  the  associated  Legendre  polynomial.  Many  papers 
have  been  published  that  tabulate  these  for  certain  values.  References  used  include 
Hobson  [27],  Snow  [28],  Siegel  [2],  Abramowitz  and  Stegun  [29],  McDonald  [30], 
Carrus  and  Treuenfels  [31],  Lebedevj  18],  and  the  Bateman  Manuscript  Project 
[17]. 

From  [17],  the  following  integral  is  the  starting  point 

nm/  /t\  1  f*  (cose +  J  Sing  cos  tr^_ 

pv  (cose)  r(l/2  _  m)  Jo  ( sint)2m  dt  (D-4) 


-dt  (DA) 


next  performing  a  change  of  variables  and  letting 


This  yields 


e}X  =  cos  6  +  j  sin  *  cos  t 


eJX 

dt  —  — ; — - dx 

sin  0  sint 


(D.  5) 


(D.6) 


and  changes  the  limits  from  0  — ►  ir  to  6  — ►  -6.  Solving  for  cos 2t,  it  can  be  found 


that  sin  t  =  Vl  -  cos2  t  and  it  is  given  by 


sin  t  —  e 


jx/2  n/2c°s  x  —  2  cos  e 
sin  6 


Putting  this  into  equation  (D.4)  with 


(cos*  +  j  sin  *  cos  t)*/+m  =  ei*(»+rn) 


(D- 7) 


(D- 8) 


and  noting  that  fg  S  =  -  f^_g,  yields 


1  (sin  fllrn  fv  eji(i/+i/2) 

P?( cos  e)  =  -=  — — r  /  - — r-rrdx  (£>.9) 

y/2n  P(l/2  -  m)  J_g  (cos  x  -  cos*)m+1/2 

This  can  be  simplified  further  since  it  is  over  a  symmetric  interval.  That  is, 
rO  cB 

f_ge}Xdx  goes  to  2/0  cosxdx  since  J_g  sin  xdx  =  0.  Therefore 
pm,. ,  _  ./j  (sin  e)m  (*  cos(t/  +  l/2)i 

K  (  1  V  »  r(l/2  -  m)  L  (COSI  - 


V  7T  r(l/2  -  m)  y0  (cos  i  -  cos*)m+1/2 
for  0  <  *  <  7r,  m  an  integer,  v  real. 


(D.10) 


BE 


which  is  Equation  (27)  in  [17].  However,  this  integral  will  have  convergence  prob¬ 
lems  for  m  >  0.  An  improper  integral  of  the  form  -£pdx  will  converge  for  p  <  1, 
and  diverge  for  p  >  1.  In  order  to  avoid  this  problem,  the  expression  is  evaluated 
at  m  =  -m  which  yields 

.  ..  [2  1  f 9  cos(i/  +  1/2)*  ,  . 

P „  (cos0)  =  \  j  —  — — — — —7 - — -  /  — — - —  7 - - dx  (£>.11) 

V  w  (sintf)mr(m  +  1/2)  Jo  (cos  x  —  cos  0)l/2~m 

but 


dx  (jD.11) 


o  coS«)  = 


(D- 12) 


So,  in  conclusion 

P?(co  s9)  = 


|2l>  +  m  +  l)  (-l)m 
V  ■nY[u  —  m  +  1)  f(m  +  1/2) (sin  9)' 


cos  [u  +  l/2)x 
(cos  x  —  cos  9)  1/2-m 


for  0  <  9  <  ?r,  m  an  integer,  v  real.  (D.13) 

which  is  the  expression  in  [18].  This  integral  has  no  convergence  problems.  When 
m  —  0,  the  exponent  is  1/2,  which  meets  the  criteria  for  convergence.  When 
m  >  0,  the  denominator  moves  into  the  numerator  and  it  is  no  longer  an  improper 
integral.  Although  the  m  —  0  case  converges,  it  is  slow.  To  avoid  this,  a  recursion 
formula  found  in  [18]  is  used  with  m  =  0.  This  yields 

P°(cos0)  =  ^  ^  2^Ft/1(cos0)  +Pj(cos<9)  (£>.14) 

Another  important  recursion  relation  is 

d  m  {v  -  m  +  l)P!J\  ,(cos0)  -  [u  -I-  1)  cos  9 P™  (cos  9) 

-prices  «)  =  i -  -  IV — - — (C15> 

which  expresses  the  derivative  of  the  Legendre  polynomial  in  terms  of  the  polyno¬ 
mials  themselves.  This  expression  is  used  in  the  determination  of  the  eigenvalues 
that  satisfy  Equation  (D.2). 

In  the  next  section,  orthogonality  is  proven,  and  the  incomplete  Legendre 
integrals  is  evaluated. 
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d  nmt  (1/  —  m  +  1)P^j(cos0)  —  (1/  +  1)  cos  9P™(cos  9) 

-P„  (cos  9)  - - - 


•V  V 

ESS® 


D.3  ORTHOGONALITY  AND  NORMALIZATION 


The  cone  presents  a  special  problem  in  terms  of  the  orthogonality  relations. 


As  seen  in  Appendix  B, 


f6°  , 

/  [-P^cos  0)]2  s\nOd0 

Jo 


(D.  16) 


The  first  thing  to  do  is  to  prove  the  original  statement  of  orthogonality  that  was 
made  in  Appendix  B.  To  this  end,  we  begin  with  the  original  differential  equation 
for  different  u  and  p .  We  will  work  with  the  form  in  terms  of  the  variable  x,  and 
then  convert  to  cos  6.  This  yields 

i  [(•  -  +  {w*  + 1)  -  ~)  JTM  - o 


Tt  [(1  -  +  {"t"  + »)  -  }  W  =  0  (O.I8) 

After  multiplying  each  equation  by  the  other  solution  and  subtracting,  what  re¬ 


mains  is 


P?(z)P?(x)  =  - 


_  d  [(»  -  *2)  ~  CTgffl)] 


W  -  p)(u  +  p  +  l) 


(D.  19) 


Next,  if  integrating  from  1  «o  ..  ^  in  x  (equivalent  to  0  to  6q  in  theta)  yields 

r>(„e-„u,  -  (,  -  10 

J 1  1  j  (*  -  mU*  +  M  +  1) 


(D.20) 

When  u  ^  p,  the  integral  goes  to  zero  at  the  lower  limit  because  of  the 
(1  -  x2)  term,  and  it  goes  to  zero  at  the  upper  limit  due  the  boundary  conditions. 
Therefore,  orthogonality  has  been  shown. 


When  v  =  fi  we  have  a  jj  condition.  Applying  L’Hospital’s  rule  in  v  yields 


/,  I^WI  ^  =  j!”!1  - *2)' 1 - iTT Ti - 


(D.21) 


Now  performing  a  change  of  variables  such  that 


i  =  cos  0 


(£>.  22) 


dx  =  -  sin  0 


(£>.  23) 


d  —Id 
dx  sin  0  rf# 


(ZJ.24) 


The  expressions  for  the  integral  are  different  depending  upon  which  eigenvalues 


are  used.  This  depends  upon  which  type  of  boundary  condition  is  satisfied.  After 


performing  the  substitutions  above,  the  integrals  have  two  forms.  When 


P?(cos0o)  =  0 


(£>.  25) 


^(cos  0O)]2  sin  0d0  =  |^j^7m(cos  *o)  ^  P™(  cos  0O) 
en 

^P™(cos0o)  =  0 

ation  integral  becomes 


(£>.26) 


(£>.27) 


os  0)] 2  sin  0d0  =  -^y  P-(cos  0O)  g~ggP^(cos  g«) 
nal  simplicity,  define 


(£>.28) 


,  define 


i™(»°)  =  £/T(cos*0) 


(£>.29) 


-.j,  JV*  l  *»  *■  f.f* 


0*0)  =  5^pr(co»*0) 


(£>.30) 


The  derivative  with  respect  to  theta  was  done  using  recursion  formulas;  how¬ 
ever,  the  derivative  with  respect  to  the  degree  was  done  using  a  double  precision 
library  subroutine.  The  expressions  are  the  same  as  those  given  in  [33].  Moreover, 
the  results  check  with  the  tables  given  in  [lj. 

In  the  next  section,  the  approach  used  to  determine  the  particular  eigenvalues 
v  and  -y  will  be  discussed. 

D.4  EIGENVALUES 


The  determination  of  the  zeros  as  shown  in  Equations  (D.l)  and  (D.2)  is 
a  very  important  aspect  of  the  solution  to  the  cone  problem.  The  solutions  to 
these  problems  determines  over  what  values  the  series  is  summed.  There  have 
been  several  articles  published  on  this  subject.  In  [31],  there  are  tables  for  both 
P™{ cos  9)  and  ^P™(cos  6).  The  results  obtained  in  this  derivation  compare  very 
well  with  those  published.  Also,  in  [2],  the  values  of  v  and  7  that  were  used  in 
their  solution  compare  very  well  with  ours.  To  find  the  zeros  begin  with  Equation 
(D.13).  The  leading  coefficients  will  not  be  zero  in  the  specified  range.  Therefore, 
finding  the  zeros  of  P™[co$0)  amounts  to  determining  where  the  integral  goes  to 
zero.  That  is, 


[**  cosfi/  +  l/2)i  . 

f(v,m,9)  =  - — — - -Ul2_mdx  =  0  (£>.31) 

Jo  (cos  x  —  cos  9y<  “  m 

as  a  function  of  1/,  for  a  given  m  and  9. 

The  determination  of  the  zeros  of  jgP™(cos9)  employs  the  recursion  formula 
in  Equation  (D.15).  The  integral  forms  from  Equation  (D.13)  were  substituted  into 
the  recursion  expression.  After  manipulating  the  integral  forms,  applying  several 
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gamma  function  identities  and  factoring  out  all  common  terms,  what  remains  is 


sp”(co5( 


1  [2  I>  +  m  + 1) 

(sin  0)m+1  V  7T  I>  —  m  +  l)r(m  H 


l)r(m  +  1/2) 


x  \(u  +  m  +  l)Ji  -  (»/  +  1)  cos  0/o]  =  0 

(D.32) 

where 

/o  =  /'  +  */*)*  ix 

Jo  (cos  X  -  COS0)1/2  m 

(£>.33) 

J  [9  cos(i/  +  3/2)a: 

1  Jo  (cos  I  -  COS0)1/2-m 

(£>.34) 

As  before,  the  leading  coefficients  are  not  zero;  therefore,  finding  the  zeros  required 

the  solution  of 

(1/  +  m  +  1)/}  —  (1/  +  1)  cos  eio  =  0 

(£>.35) 

as  a  function  of  t/,  for  a  given  m  and  0.  These  results  also  compared  very  well  with 
the  published  values.  All  of  the  preceeding  integrals  were  performed  using  a  32 
point  Gaussian  quadrature  routine.  The  roots  were  found  using  existing  library 
subroutines. 

In  the  last  section,  the  method  employed  to  verify  the  validity  of  the  Legendre 
function  routines  is  discussed. 


D.5  VERIFICATION 


In  [33],  a  Wronskian  result  is  given.  A  more  general  form  of  this  expression  is 

P^{cosei)8in02^P^{cos92)  -  P^cos^Jsintf^P^cosfli) 

=  }  sin[(i/  -  (D.36) 


Letting  02  =  x  -  9\y  then 


d0  2 


P?(cos02)  =  -^-P™(~cos0,) 


(£>.37) 


Therefore,  Equation  (D.36)  becomes 


—  _  2  sin|(t/-m)irl  r(t/+m+l) 

—  »r  8in0j  r(t/-m+l) 


(D.38) 


The  routines  were  checked  against  existing  tabulated  results;  however,  for 
angles  and  degrees  that  were  not  tabulated,  Equation  (D.38)  was  used  to  verify 
the  routines.  This,  in  effect,  checked  both  the  routine  for  P™(c os  6)  and  the  routine 
for  4$  P™ {cos  6)  simultaneously. 
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